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Abstract 

We discuss within the weak-field approximation and the derivative expan- 
sion how the area law of the Wilson loop follows directly from the vacuum 
condensate of mass dimension 2, i.e., simultaneous Bose-Einstein condensation 
of gluon pair and ghost-antighost pair. Such a novel vacuum condensate was 
recently claimed to exist as the non-vanishing vacuum expectation value of a 
BRST-invariant composite operator of mass dimension 2. First of all, we use 
a version of the non-Abelian Stokes theorem to rewrite the Wilson loop line 
integral to a surface integral. Then we convert the Yang-Mills theory with an 
insertion of the Wilson loop operator into a bosonic string theory with a rigidity 
term by way of an equivalent antisymmetric tensor gauge theory which couples 
to the surface spanned by the Wilson loop. This result suggests an intimate 
relationship between quark confinement and mass gap in Yang-Mills theory. In 
fact, the dual Ginzburg-Landau theory describing the dual superconductivity 
is also derivable by making use of duality transformations without using the 
naive Abelian projection and without breaking the global color invariance of 
the original Yang- Mills theory. This feature is desirable from the viewpoint of 
color confinement preserving color symmetry. 
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1 Introduction 



Proving quark confinement is still one of the most important problems in theoretical 
physics. Based on intensive investigations in the last decade, the dual superconductor 
picture for the vacuum of quantum chromodynamics (QCD) is believed to be the 
most promising scenario of deriving quark confinement. To derive this picture, it is 
convenient to select out an Abelian sector of the original Yang-Mills theory so that 
the duality transformation of the Abelian sector gives a dual Abelian gauge theory 
in which the dual superconductivity is caused by magnetic monopole condensation. 
Such a procedure is called the Abelian projection 0. Therefore, the Abelian projected 
effective gauge theory of the original Yang-Mills theory has been intensively studied in 
the last decade from the viewpoint of numerical |3|-|7| and analytical |8HlO| methods. 

Nevertheless, we wish to point out in this paper that the naive procedure of the 
Abelian projection is not necessarily indispensable to derive the dual superconduc- 
tivity for explaining quark confinement. Even if we begin with the Lorentz gauge 
dfx^ix) = which has global gauge invariance, the dual superconductivity can be 
derived without breaking the original global gauge symmetry, as will be demonstrated 
in this paper. This result should be compared with the Maximal Abelian (MA) gauge 
which breaks the global gauge symmetry from the original gauge group to the maxi- 
mal torus group (except for the discrete Weyl symmetry JTTJ ) . This choice of Lorentz 



gauge can greatly simplify the actual calculations by virtue of preserving the global 
gauge invariance. This fact is also desired for enabling us to tackle color confinement 
problem (i.e., confinement of all color non-singlet objects) beyond quark confinement 
problem. We want to add that we do not deny the MA gauge as a choice of gauge 
for deriving quark confinement, since quark confinement should be a gauge invariant 
phenomenon. Thus we conclude that the dual superconductivity of QCD vacuum is 
an intrinsic property, independent of the choice of the Abelian projection to define 
the monopoles. 

In this paper we adopt the Becchi-Rouet-Stora-Tyutin |L2] (BRST) invariance as 
the first principle to characterize the quantized Yang-Mills theory, rather than the 
gauge invariance which is broken by the procedure of gauge fixing in the course of 
quantization. The main purpose of this paper is to demonstrate how to derive a con- 
fining bosonic string representation of Yang-Mills theory. In this derivation, we show 
that the Wilson loop average exhibits the area law decay if the vacuum condensation 
of mass dimension 2 occurs in Yang-Mills theory, i.e., the vacuum expectation value 
(VEV) of the composite operator [ |13|j , 

1 

2" 



o = n- 1 J d 4 x tr G/H 



-af"{x)^{x) + \i c i{x) c i{x) 



XD 



is non-vanishing (O)ym 7^ 0, where Q is the volume of the space-time Q :— J d 4 x and 
A is a gauge fixing parameter in the Lorentz gauge d^^/^x) = 0. 

The vacuum condensates of mass dimension 2 were recently proposed by several 



authors [|T4] -|T7H and the physical implications have been extensively studied in a 



couple of years [[T3|, P~8| pOfl - Recently, it has been shown JT3[ that the composite 
operator O of mass dimension 2 can be made both BRST and anti-BRST invariant^] 



1 Here the BRST should be understood as the on-shell version (2.5) of the BRST transformation 



which is obtained by eliminating the (Nakanishi-Lautrup) auxiliary field 38 from the off-shell BRST 



transformation (2.4). Note that the operator O does not include the SS field. 
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in the generalized Lorentz gauge ET] and the modified MA gauge ||22[| . Especially. 



in the limit of Landau gauge A — > 0, the vacuum expectation value of O reduces to 



the gluon condensation of mass dimension 2 proposed in 23 , although the way of 



taking the limit is not unique, see pi . Therefore, we claim in this paper that quark 



confinement follows from the vacuum condensate of mass dimension 2'm. a BRST and 



anti-BRST invariant manner, as suggested in the previous paper |fL3 |. 

In this paper, we emphasize the importance of a novel gluon condensate of mass 
dimension 2, in contrast with the conventional gluon condensate of mass dimension 
4, i.e., {J-? u ) which is obviously gauge invariant [25| and BRST invariant.f] In fact, 



our results indicate an intimate connection between the existence of mass gap and 
quark confinement in Yang-Mills theory, since the non-vanishing vacuum condensate 
of mass dimension 2 leads to the dynamical generation of the effective masses of 
gluon and ghost. In order to elucidate this statement, we rewrite the Yang-Mills 
theory in the presence of a Wilson loop operator into a bosonic string theory within 
the weak-field approximation. It turns out that the string action is given by a sum 
of the Nambu-Goto term and the rigidity term (as announced in p7||). Therefore, 



the resulting string is called the rigid string hereafter, following Polyakov |36| . |37 
also PP. 



sec 



We do not attempt to analyze in detail the physics of the obtained rigid string, 
since such investigations have been performed by numerous authors in recent years. 
See e.g., P^-BSl for the old works of string representation of QCD and p]|-^T| for re- 



cent developments on this subject. The novel result of this paper is the determination 
of the parameters in the rigid string action in terms of the original Yang-Mills theory 
(within the approximations adopted). Therefore, their result can also be applied to 
our string theory by substituting the determined parameters in our paper into their 
results without repeating the calculations. The main purpose of this work is to build 
bridges between Yang-Mills theory and a rigid string theory based on the existence 
of a novel vacuum condensate. We show that an antisymmetric tensor gauge the- 
ory (ASTGT) of rank 2 plays the dominant role in establishing the relationship, as 
expected in the previous investigations. 

Our analysis begins with a recent observation that non-perturbative corrections to 
perturbative results do indeed exist even in the high-energy region. Such corrections 
can lead for example to the linear static potential for a pair of quark and anti-quark 
at short distance, although the linear static potential is expected to occur in the 
long distance according to the conventional wisdom.^ Consequently, a short confining 
string will be obtained as the string representation of the Yang-Mills theory in the 
relevant energy region (except for the extremely low-energy). Extremely low-energy 
region is still beyond our approach due to lack of suitable tools of analysis. We hope 
that the strategy given in this paper will survive in the low-energy region, if the 
suitable tools for the calculation become available. 

This paper aims to demonstrate that quark confinement can be understood at 
least qualitatively even in the lower order of approximations (expansions) adopted, 
just as the area decay law of the Wilson loop average can be derived even in the 



2 The derivation of a confining string theory due to gluon condensate of mass dimension 4 has 
been already tried by several authors, see e.g. a review by Antonov p(| and references therein. 

3 Therefore, our study on quark confinement does not imply the existence of 1/p 4 infrared sin- 
gularity in the gluon propagator which was investigated e.g. by Ellwanger pa]. 
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lowest order of the strong coupling expansion in lattice gauge theory. In lattice gauge 
theory, it is necessary to control the continuum limit for the complete proof of quark 
confinement. In the continuum formulation, this problem is absent, but we must 
overcome other difficult problems to remove the approximations. 

The outline of this paper is as follows. In section 2, we enumerate the necessary 
steps to arrive at the main results. In section 3, we give some remarks on the results 
obtained in section 2. In the final section, we give conclusion and discussion. In order 
to clarify the essence of this work, technical details are all omitted in section 2. They 
are given in Appendices A to F. 



2 Steps to arrive at the main results 

The main results of this paper are obtained following the steps enumerated below. 

2.1 Step 1: Definition of the Wilson loop operator and its 
average 



The Wilson loop operator Wc\&* \ is defined by the trace of the path-ordered product 
of the line integral of the non-Abelian gauge field ^ i.e., Lie-algebra Q valued one- 
form = stf^dx^ = srf^T dx^, along a closed loop C: 

WcW] = A^trjpexp ig j> dx^stt^x) j, (2.1) 

where M is the normalization factor to guarantee W^clP] = 1 and we adopt the con- 
vention in which the coupling constant g is explicitly written (For a precise definition 
of the path ordering V, see e.g., [0]). In what follows, we assume that the loop C is 
a connected closed path without self-intersections. 

The Wilson loop average W(C) is defined as the vacuum expectation value (VEV) 
of the Wilson loop operator Wc\^\ in the Yang- Mills theory. For concreteness, we 
suppose that the manifestly Lorentz covariant formulation of the Yang-Mills theory is 
adopted. In this paper, we always work in the Euclidean formulation, unless otherwise 
stated. In the functional integral formulation, the Wilson loop average is defined by 

W(C) = {WcW\)ym = Zylt J dfi YM e- s ™W c [^], (2.2) 
where dfiyM is the measure of functional integration, 

dfi YM = VsrfpStfDWtf (2.3) 

defined as the product measure of the respective field, i.e., the gauge (gluon) field 
sd^ Nakanishi-Lautrup (NL) auxiliary field Faddeev-Popov (FP) ghost field ^ 
and antighost field ^ ', and Sp M is the total action of the Yang-Mills theory consisting 
of the pure Yang-Mills term Sym for a gauge group G, the gauge-fixing (GF) term 
Sqf and the associated FP ghost term Spp, i.e., SpM — Sym + Sqf + Sfp- The 
gauge- fixing part Sgf+fp '■= Sqf + Sfp is determined based on the Becchi-Rouet- 
Stora-Tyutin [12|] (BRST) transformation 6b- The BRST transformation in Euclidean 
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space is given by 



S B ^{x) = f^KMx) := d^(x)+g(^(x) x <#(x)), 
1 



Sb^(x) = -I 
8 B 38(x) = 0. 



-g(?(x) x tf(x)), 



[x 



(2.4a) 

(2.4b) 

(2.4c) 
(2.4d) 



In fact, the gauge- fixing part Sgf+fp '■= Sgf + Sfp is written in the BRST exact 
form p5| , i.e., Sgf+fp = ^b( - • • ) using a quantity (• ■ ■ ) depending on the gauge-fixing 
condition. Therefore, the explicit form of Sgf+fp depends on the gauge fixing con- 
dition. The pure Yang-Mills action is BRST invariant, SbSym = 0. The gauge-fixing 
part is BRST invariant, SbSgf+fp = due to nilpotency of the BRST transformation 
5| = 0. Consequently, the total Yang-Mills action is BRST invariant, SbSym = 0- 

For our purposes, however, it is more convenient to eliminate the NL field 38 by 
integrating out it in the beginning, since the composite operator O defined by ( |1.1|) 
does not depend on 38. Even after the elimination of 38, the total action is invariant 
under another BRST transformation which we call the on-shellBRST transformation. 
The explicit form of the on-shell BRST transformation depends on the choice of the 
gauge fixing condition. In the generalized Lorentz gauge |23j, the on-shell BRST 
transformation reads 



Sb^(x) 
Sb^(x) 



- l -g{^{x) x V{x)), 



~d^(x)+Zgtf(x)xtf(x) 



(2.5a) 
(2.5b) 

(2.5c) 



where A and £ are two gauge fixing parameters in the generalized Lorentz gauge. 
In the £ = 0, the generalized Lorentz gauge reduces to the ordinary Lorentz gauge 
as found in the textbook of quantum field theory where a gauge fixing parameter 
A = corresponds to the Landau gauge. It has been shown [|13| that the composite 
operator O defined by ( |1 . 1|) is on-shell BRST invariant, if the theory is restricted to 
the subspace specified by A = or £ = 1/2. It is confirmed that this property is 
preserved under the renormalization and that the renormalization group flow is 
restricted to the subspace if it starts at a point in the subspace (In particular, A = 
is a infrared fixed point). In what follows, we restrict our analysis to the first case of 
A — > with £ = for simplicity. The second case will be discussed elsewhere. 



2.2 Step 2: Non-Abelian Stokes theorem 

First of all, we attempt to rewrite the non-Abelian Wilson loop operator into a quan- 
tity written in terms of a surface integral over an arbitrary surface S with a boundary 
C, i.e., <9£ = C (See Fig. [[]), by making use of a version []64j of the non-Abelian Stokes 
theorem (NAST) ||65|| . In this paper, we adopt the Diakonov-Petrov (DP) version [64 



of NAST which does not contain the ordering of the surface integral, although there 
are a number of versions of NAST, see references cited in ||62|| . The DP version of 
the NAST was first obtained for G = SU(2) and the same result was rederived for 
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Figure 1: Wilson loop C and the surface E. 



G = SU (2) by making use of the coherent state representation |)6| and later extended 
to SU(N) by making use of the generalized coherent state representation |S2|,|67|. For 
the simplest case of G = SU(2), the NAST is given by 



W, 



c 



dns{V) exp 



J 



v fV 



(2.6) 



where a character J specifies the (2 J + l)-dimensional representation to which the 
probe source of quark for defining the Wilson loop operator belongs (J = 1/2 for 
the fundamental representation, J = 1 for the adjoint representation, and so on). 
Here dusiV) is the product measure Yixes dfi(V(x)) of the invariant Haar measure 
d^(V(x)) on the coset G/H = SU(2)/U(1)(= S 2 3 V) which is defined at an arbitrary 
point x on the surface S. For an infinitesimal surface element dS^, an antisymmetric 
tensor field f^ v is defined byf] 



f, 



X 



9 M [V(x) • K{x)\ - d v [V{x) ■ ^{x)\ - g- x V(x) ■ (^V(x) x d v V(x)) 



(2.7) 



where V is the flag variable V(x) := V (x)a /2 with Pauli matrices a (A = 1, 2, 3) 
and three components of V^ A (A = 1, 2, 3) constitute a unit vector, i.e., V(x) • V(x) = 
V (a;)V" (a;) = 1. Here the last term in (|2.7|) corresponds to the topological term. 
By parameterizing a unit vector V explicitly as V\ = sin 9 cosy?, V 2 = sin^sin^, V3 = 
cos 8, the invariant Haar measure on SU(2)/U(1) reads 



d/x(^(x)) = 6(V(x) ■ V{x) - l)dV 1 (x)dV 2 {x)dV 3 (x) = sin 6 (x)d6(x)d<p(x). (2.8) 



eV 

1 flU 

for Surtax 



The /Y, is invariant under the non-Abelian full gauge transformation, i.e., 5^,/ 



V fl uj(x) and <5 w V(x) 



x 



ig[u;(x), V(x)] in the infinitesimal form. 



4 This tensor has the same form as the 't Hooft tensor for describing the 't Hooft-Polyakov 
magnetic monopole, ifV A (x) is identified with a unit isovector of the Higgs scalar field, (f> A :— (j) A /\(f>\. 
This fact suggests a possible interplay between magnetic monopole and quark confinement. In other 
words, the magnetic monopole might be encoded in the Wilson loop operator, see [E3,p§,pq,Q. 
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It is invariant also for a finite gauge transformation: ^(x) — > f/(x)^(a;)L r ^(x) + 
^-^^^^(s), and V(z) -> C/(a;)V(x)t^(x) for C/(x) = e iflW ^ with = 
ui (x)cr /2. Similarly, the BRST transformations of V is defined by replacing u 
with the ghost field i.e., 

S B V(x) = ig[tf(x), V(x)] = gV(x) x <g(x), (2.9) 

just as the BRST transformation of ^ is given by Sb^(x) = &^€{x). Then the 
Wilson loop average reads 



W{C) = Zy l M J dfXYMe'^M J dfi s (V)ex P [z<^ JjS^f t 



[2.10) 



2.3 Step 3: Introducing auxiliary antisymmetric tensor field 

Next, we introduce an antisymmetric tensor fieldf] h^ v (x) (of canonical mass dimension 
2) by inserting an identity into the functional integral: 



J exp { - JL J d 4 x[h^(x) - /^(x)] 2 } , (2.11) 



where k q is an arbitrary dimensionless constant and the normalization factor is in- 
cluded in the measure^ T>h^ u . The simplest way to keep the gauge and/or BRST 
invariances is to require the relation, 5 u [h^ v {x) — fY v (x)\ = and/or 8-Q[h^ u {x) — 
fuv( x )\ = 0- Hence, the BRST invariance is preserved, if h^ v is an Abelian tensor 
invariant under the BRST transformation, i.e., 

«bMi) = ) (2.12) 

since S B f^(x) = 0. 

It should be remarked that the tensor field h^ v is an auxiliary field without its 
kinetic term at this state. However, the kinetic term will be generated by radiative 
corrections in the quantum theory. This is one of the most important points in 
deriving the string representation of the Yang-Mills theory, as will be discussed in the 
following. 

For ( [2. 11| ) to be meaningful, furthermore, V(x) must be defined on the whole 



space R 4 . Therefore, the field \[x) on the surface £ must be continued to the 
whole space R 4 outside the surface S. This continuation would not cause a trouble, 
since this continuation does not lead to the failure of uniqueness due to encountering 



5 The antisymmetric tensor h^ v has no index for internal symmetry. It is possible to introduce 
the antisymmetric tensor field h^ v in the adjoint representation of the gauge group G. However, 
we do not use this type of tensor. The reason is as follows. The tensor is not gauge invariant. 
Therefore, the Abelian components are identified with the diagonal components and selected out by 
making use of a naive procedure of Abelian projection, see Ellwanger ]5j|. However, this procedure 
breaks the global color symmetry. 

6 For the integration by parts to be possible, the translational invariance of the measure is assumed 
to hold, i.e., V(h + a) = T>h, which implies an identity, J Vhj^-f(h) — for arbitrary functional 
f(h) of h. This is in fact necessary to define the Gaussian integration Q2.ll ). The constant kq will 



be determined so as to reproduce the correct normalization of the dual Ginzburg-Landau theory in 
the final stage of the derivation. 
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the singularity (Such a possibility was carefully examined when the line integral is 
continued into the surface integral in NAST by Diakonov and Petrov [fTOD J 7 ] Therefore 
we can replace dns(V) with dfj,jn(V) := dfi^V) without any difficulty (This issue 
will be further discussed in a forthcoming paper). 
Thus the Wilson loop average is written as 



W(C) = J Vh^exp^— J d A xhl v (x)}Zy l M j dp, 



x 



exp [1| d'x - lCC)}<*P [i9^ f^f, 



otot 

YM e a ™ 



J 



rV 

' jJLV 



where we have introduced the expectation functional, 



(2.13) 



([■••]>v:= d^{V)[ ■ ■]■ 



(2.14) 



2.4 Step 4: Cumulant expansions 

By introducing the vorticity tensor current Q fJlV (x) by 



d 2 S llu (x(a))S 4 (x - x(cr)), 



(2.15) 



it is possible to cast the surface integral into the space-time integral according to 
J s dS ,J ' u (x(a))f}[ u (x(a)) = J d 4 xQ^ u (x)f^ u (x), where Q^ u (x) has its support on the 
two-dimensional surface £ (parameterized by the coordinate a = (01,02)) with a 
target space coordinate x M = x M (er). We attempt to convert the interaction term of 
the surface element do~^ v with the H v into that of the surface element do~^ u with h 
In fact, an identity: 



exp 
exp 
exp 



J 



ig 



[ d 2 S^fl 
J s 

igJK I d A xQ 



exp 



1 

K 



eV 



exp 



1 

K 



_5_ 

ilu 6h 



fJ,V 



J d x ^ hfj,i/fi 
J d 4 x^h flu f t 
— f d 4 x-h UL/ f, 

KnJ 2 



(2.16) 



and the integration by parts allows us to achieve the desired result: 



W(C) 



Vh 



[IV 



X %YM 



exp 



igJn J d 4 xQ 



djJ,YM e 



exp 




(2.17) 



It should be remarked that the V field is an important but artificial field (without 
any direct physical relevance) which has been introduced in the NAST to rewrite the 

7 The surface integral is obtained from the line integral of the diagonal component := tr[(73 1 e^] 
of s&u := U^W + g-HUdpUl by the usual Stokes theorem as ^dx^a^ = J^dS 1 * 1 ' fj v . Hence 
V A {x) = tr[W(x)a 3 U{x)a A ] for U[x) = exp(ix(x)a 3 /2) exp(i6(x)<T 2 /2) exp(iip(x)a 3 /2) G SU(2), 
where Xi 9, <P are Euler angles. Note that V A is invariant under the residual U(l) gauge transforma- 
tion and hence it is the flag variable [p2[ belonging to SU(2)/U(1). See [p2,pffl for details. 
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line integral into the surface integral. Now we attempt to eliminate (or integrate out) 
the V field to obtain an effective theory written in terms of the auxiliary field h^ v 
and the fundamental field $ = {^, ^, ^}. In the formal level, this step can be 
realized as 



W{C) =Z Y l M J djiYM J Vh^exp igJn, J d 4 xQ 



5 



X 



6XP {"4^ / dAxk '^ ~ ^™ [$] + F[K ^ ] ) 



h'=h 



(2.18) 



where we have introduced F[h, sf\ which is written as the formal power series in k 1 
by making use of the cumulant expansion for the expectation value ((■ ■ • ))y: 



F[K *\ = In ( exp {1 / d 4 x (^h^f* - \f^f t 



v t v 



oo -n 

= E 



n=l 



d A x-h f v 



v 

I i n \ con 

d 4 X - fV f v 



(2.19) 



with the connected expectation value ((•••))v 11 defined from ((•••))v m the usual 
way. 

Subsequently, integrating out all the fundamental fields <£> leaving h untouched, 
we arrive at the effective theory written in terms of an antisymmetric tensor field h^ u 
alone. In addition to the first cumulant expansion for ((• • • ))y: 



W{C) = J DV (exp igJn J d A xB 



_5_ 

llu Sh 



YM 



(2.20) 



the second cumulant expansion for ((• • • ))ym results in 



W(C) = J £>Vexp igJKoJ d A x& 



S 



-s h 



h'=h 



S h :-- 



4k 



d xh 



OO 1 



171=1 



ml 



YM 



(2.21) 
(2.22) 



where ((• • • ))ym denotes the connected expectation value in the Yang-Mills theory. 

For the purpose of obtaining the effective theory written in terms of h alone 
without breaking the BRST invariance, it is better to use the form: 



W{C) = J £>V (exp igJn J d 4 xB 



_6_ 

au Sh 



\1V 



exp 



X 



V fV 



VI YM 



(2.23) 



Therefore the cumulant expansion reads 



W{C) = J PVexp igJ Ko J d A xQ, v ^- exp^-^-J d 4 xh% u + Q[h]} 



h'=h 

(2.24) 
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by making use of the cumulant expansion for the double expectation (((■ ■■))v)ym 
defined by the formal power series in Kq l : 



v t v 

jlV 



G[h] =ln (( exp {1 J d*x (± h^f* ~ \fU 
J d x—h^ u f^ v — J d x—f fMU f f 



E 

n=l 



V f V 



VI YM 
n \ com con 



V I YM 



(2.25) 



Furthermore, it is possible to calculate the derivative with respect to h in the 
closed form,F| 



exp 
: exp 



igJn / d 4 yQ^(y) 



4k 



d 4 xh 2 (x) - -g.J / d 4 xQ l _ tu (x)h IMU (x) + — (gJ) 2 no / d 4 xQ 2 (x) 



24 



(2.27) 



This procedure is not necessarily an indispensable step, but it is useful to avoid the 
complicated calculations in the later stage of integrating out all the fields except for 
h. Thus we obtain 



W(C) = J Vh^e' Sh+ ^ {9j)2 1 d4xe i^ x) exp -ig- J d 4 xG^ u h 



(2.28) 
(2.29) 



In this step, we are interested in only the terms coupled to the h^ u field, since other 
terms merely give /i^-independent constants after taking the double expectations or 
double cumulants, ((• • - ))y n and ((• • -))ym- The full form of the action written in 
terms of $ and h is quite complicated but necessary for quantitative analysis and will 
be given in a subsequent paper. 

Retaining only the terms up to quadratic in h in the cumulant expansion (which 
we call the weak field approximation |36| or the bilocal approximation)^ we obtain an 
expression of the effective theory written in terms of the h alone (see Appendix [A] for 

8 We have used the Baker-Campbell-Hausdorff formula 



e x e Y 



{X + Y + l[X, Y] + [X, Y]] + [Y, [Y, X}}) 



(2.26) 



9 It turns out that the quadratic terms are sufficient to reproduce the London limit in the type 
II superconductor. However, we can show that the quartic term in h is necessary to go beyond the 
London limit and to reach the border between type II and type I from the type II, see Appendix 
Especially, the truncation of the series excludes the possibility for the series in h to be resumed 
into a closed form which might be a periodic function of h. If so, it will be possible to distinguish 
the half-integer J representation from the integer J one in the calculation of the string tension, as 
realized in compact QED |3(|. 
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details of calculations): 

S h = [ d 4 x—h lxu (x)h fll/ (x) - —z [ (fx [ d i yh llu (x)h pa (y)g' 2 (Q pu (x)Q pa (y)) 



con 
V 



1 



d A xJ A(^ A (x)V B (y))r(<(^)< B (y))y A n /^^ a (x)^^ /3 ( Z/ ) + • • • , 

(2.30) 

where we have introduced an antisymmetric tensor of canonical mass dimension 2: 
JVC*) : = V 0) ■ (Wz) x 9 v y(x)) = -tt ni {x). (2.31) 

2.5 Step 5: Gluon condensation and dynamical generation 
of the kinetic term for the antisymmetric tensor field 

First, we assume the translational invariance for the VEV ((• • • ))y. Hence we have 

(V A (x)V B (y))r = G AB (x - y) (2.32) 

with a normalization G AA (0) = 1. In addition, we can put 

(n^(x)n pa (y))<$» = K l j tl/pa (x - y), (2.33) 

where K^ vp(J has the symmetry: K pupa = -K uppa = -K pvap = K papv . Note that 
(V A (x)) c J n = 0, and (n pu (x)) c J n = 0. 

Second, we expand h(x) and h(y) in powers of the relative coordinate r M := (x—y) p 
around the center-of-mass coordinate X := (x+y)/2. Since x = X+r/2, y = X — r/2, 
we have 

{ M$} = ^" (X) ± if + ^A^VW + ■ ■ • • ( 2 - 34 ) 

Third, we take into account the decomposition of the double integration over the 
space-time: 

J d 4 x J d 4 y{- ■■) = J d A r J d 4 X(- ■ ■ ), (2.35) 

where the Jacobian for change of variables is equal to one. 

The first term of Sh ( |2.30| ) is the mass term for h pv . The second term of Sh ( |2.30D 



is cast into 

J d A x J d 4 yh pu (x)h pa (y)(Q pu (x)Q pa (y))^ 
= J d A rK pupa {r) J d A Xh pv {X)h pu {X) 



+ J d A rK pvprj {r) l -r ar!i J d A Xh pu {X)d a dph p(J {X) + 0(hd 4 h) 

+ Ifa® J d A Xh pu (X)I pup(T d 2 h p(T (X) + 0(hd A h), (2.36) 
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where we have put 



d\{r 2 ) l K t 



pupa \ 



fa® I 



pupa j 



L pupa 



3 pa & up) 



(2.37) 



In the tree level (without loop corrections), the kinetic term for the field V A (x) is ab- 
sent.^ Hence V (x) is subject to the short-range correlation,^ (V A {x)V B {y)) c J n = 
CS AB S 4 (r), where C must be a constant with mass dimension —4 based on dimen- 
sional analysis, i.e., C = C Aj^ with a dimensionless constant C (In other words, 



/ d A x is the volume of space-time). The correlation of 



con 
V 



c't s 

^ ± pupa u 



is also ultra-local: 
where C is another dimensionless constant. There- 
for £ > 1 in the tree level. This implies that the 



(tt pu (x)tt pa (y)) 
fore, we have a/ -* ^ and 
kinetic term for h is not generated from the second term in the RHS of Q2.30 ) in the 
tree level. The contribution from renormalizes the mass term and can be 

absorbed by the redefinition of Kq- 

It is possible to relate the action ( 2.30|) to the conventional gluon condensate of 
mass dimension 4 [p^[] , as will be discussed in the next section. On the contrary, we 
discuss in this section how to relate the action (|2.30|) to a novel vacuum condensate of 
mass dimension 2 



expansion (OPE):Q 



17]. A pair of gluon field operators has the operator product 



)i + vy [ 5 



AB i 
pu \ 



}AB i 



1 

r 

+ 



(2.38) 



where the first term D AB is the gluon propagator dressed by the perturbative quan- 



tum corrections with the tree expression given by (D ) AB (r) = 5 AB (l/d 2 )[5 pu — (1 — 
X)d p d u /d 2 ]5 A (r). 

In this section, we focus on the Lorentz gauge condition d^stf^x) = as a mani- 
festly Lorentz covariant gauge fixing. (Other types of gauge fixing will be examined 
in the forthcoming paper.) In the Landau gauge A = 0, the Wilson coefficient has 

W [^^]AB Q) whne W W-g\AB = qJT3| Thig result wag 



been calculated in 16, 19 



also confirmed as a very special limit in the framework of the most general Lorentz 



10 In order to make the expectation value ((• • • ))y well-defined, however, we need the kinetic term 
of the V field which comes out in the course of the derivation of NAST and goes to zero in the limit 
of removing the regularization, see for details. Such a kinetic term can also be generated by 

radiative corrections, as will be examined in a subsequent paper. 

11 The support of the relevant correlation function consists of a point x = y, if the kinetic term is 
absent. Therefore, it must be proportional to the delta function (and its derivatives). This is because 
for x ^ y, the absence of the kinetic term implies (V A (x) V B (y))^ n = (V A {x)) c J n (V B (y))^ on = 0, 
and {n^(x)n pa (y)) c v ™ = (n^(x)} c v ™(n pa (y)} c v ™ = 0. 

12 It should be remarked that the gluon condensate of mass dimension 4, i.e., {^ v (x)) does not 
appear in the OPE of the gluon propagator, see Q). Therefore, we cannot use it as a basic ingredient 
to derive the mass gap and quark confinement in our strategy. 

1 In gauges other than the Landau, the operator mixing among the operators with the same 
mass dimensions and the same symmetries does take place so that W ^ 0, as demonstrated 
in H). Only in the Landau gauge A = 0, the operator mixing disappears. The ghost-antighost 
composite operator mixes with the gluon pair composite operator under the renormalization. Here 
a BRST invariant combination of two composite operators appears in the calculation beyond the 
tree level. See Hi for details. 
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gauge fixing |13],|23| where W^jjf ^ in general. For SU(N C ) Yang-Mills the- 
ory, the RG improved value of the Wilson coefficient is given by W^^ip) = 
S AB T^ U (p) ( ) 4 ^7 j n momentum representation where /?o := y-^c and T^*" is the 
tree result: Ti V {p) = 2 (jv| 9 -i) [(A)) 2 -Pr]^ with the transverse projection P^, := 5^ — 
PiiPv/p 1 - Here A is the RG invariant constant in Yang-Mills theory corresponding to 
K QCD in QCD. Therefore, W^^ip) = < J^_^L_^ ( ^ ) ! fpT _ Qn the other 



hand, the RG improved gluon propagator is given by D^{jp) = ( |~^ ) 6 ^ (D ) AB (p), 



sec 



19,24 



Hence the third term of ( p. 30 ) reads 



rf 4 x y d 4 y(y A (x)^(y))r(<(^)< < (y))r M^V(^)^^(2/)> 
d 4 x d a h m {x)—dph ll p(x) 

1 OO „ ,1 > 

- -J2 K $j ^ 4 X(-<(X))^9 Q ^(X)(9 2 )^^(X), (2.39) 

where we have defined 

K$ = a e [ d\ (r 2 YG AB (r)W^ A J(r), (2.40) 



with appropriate numerical coefficients ag{£ = 0, 1,2, •••), in particular, a = 1. 
Note that the canonical mass dimension of Kf® is — (4 + 2£). The first term of the 
right-hand side (RHS) of ( |2.39|) comes from the first term of RHS of ( 2.38Q (without 



perturbative radiative corrections) and contributes to the renormalization of the mass 
term, leading to a modification of Kq. Hereafter the vacuum condensate {\&/p(X)) 



YM 



in ( 2.39[) should be understood to be an Euclidean quantity together with the Wilson 



coefficient W$**lj£ in ( ggOp 



The first two terms £ = 0, 1 in the series are sufficient for our purpose at moment. 
Note that Kj$ ^ and Kj® = for £ > 1 in the tree level. The term proportional to 
c</°) is absorbed into the redefinition of Kq. Therefore, we have only to consider Kjff). 
Thus we obtain the effective theory (up to higher power terms and higher derivative 
terms) : 

s h = [ d 4 x-^/v(x)/v,(x) - -L ( d A x d a h m (x)^d p h^(x) 

+ ^E^ / d'xl\^{X)\ d a h, a (X)(dyd (3 h u(S (X) + ... . (2.41) 

Even in the tree level, therefore, the kinetic term for h is generated from the third term 
of (p.30Q if the vacuum condensate of mass dimension 2 exists as a non-perturbative 
effect .0 Thus we have observed that the kinetic term of h^ v is generated by the 
existence of vacuum condensate of mass dimension 2. 

14 On the other hand, if V A (x) was subject to the short-range correlation as in the tree level, the 
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This phenomena was already pointed out to occur in the Maximal Abelian (MA) 
gauge in the previous paper |HJ where the ghost-antighost condensation |fH|, |1"5| 
played the similar role to the gluon pair condensation. In the previous paper, however, 
we have used from the beginning the massive propagator for the off-diagonal gluon in 
the MA gauge in the course of loop calculation. In such a case, quadratic divergence 
may appear. The derivation just given in this paper does not use such an assumption, 
although it is expected that the gluon pair condensation yields the massive gluon. 
Therefore, we can avoid the problem of quadratic divergence in our approach. In 
fact, we can explicitly show that there is no quadratic divergence due to cancellation 
among Feynman diagrams (Fig.l, Fig.2) indicated in @, provided that the BRST 
symmetry is not broken. 



2.6 Step 6: Effective antisymmetric tensor theory 

In order to obtain a dual theory of the Yang-Mills theory, we introduce the dual 
variable B pv as the Hodge dual of h pu , i.e., 



B 



*h = _ e hp° 



It is obvious that B^ u is both gauge and BRST invariant. Note that 

^^pdah^adphyp =-(e pupa d u B pa ) 2 . 



(2.43) 

(2.44) 
(2.45) 



Similarly, we have for any integer 



■1,0,1,... 



6^d a h m (cP)%Kp = -(e^Bn^Yie^^B^). 



(2.46) 



Thus we arrive at a BRST-invariant dual theory which is equivalent to Yang-Mills 
theory with an insertion of the Wilson loop operator: 



W(G) ^ZyltJvB^e-^+^^f^U^exp -ig J - J d'x *Q pu B 



(2.47) 



S d [B] 



1 r>2 , 

4k » v 



d 4 x 

g4 R 2 BfivI/xvpaO B pa 



^{f»e°d v B pa f - ^^ vpa d v B^)^{e m ^B^) 



32 K/Q 



\d a e^°d v B pa f + 



(2.48) 



effective theory in h was simply obtained without relying on the OPE as 



d^x 



1 h^ l/ (x)h l _ iy (x) + ^^^.s/ p (x))^^d a h pa (x)dfjh u p(x) 



4k 



Kq 4 \ 2 



YM 



(2.42) 



where C must be a constant with mass dimension —4 based on dimensional analysis, i.e., C = CqA^ 



with a dimensionless constant Cq and a dimcnsionful scale A/^. This case is derived in (2.30) by 

)v n = C Ajfr 



taking the correlation: {V A {x)V B {y)) c ^ n = C Aj^S AB 6 4 (r). 
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where we have put 



(2.49) 



by introducing dimensionless constants K™' of canonical mass dimension — (4 + 2£), 
and the vacuum condensate a vc of mass dimension 2 is defined by the VEV of a 
composite operator O: 



(Jr. 



(O) 



YM- 



(2.50) 



For the conventional Lorentz gauge, the vacuum condensate has been evaluated only 
in the Landau gauge, see [72] for the analytical result and [H], 20] for the numerical 
results. Both results conclude that the non-zero vacuum condensate a vc does indeed 
exist where a vc = (|^ 2 (0)) ^ for the translational invariant vacuum. It is possible 
to calculate the vacuum condensate in the most general Lorentz gauge, see [O. 24 . 



2.7 Step 7: Area law of the Wilson loop average 



We proceed to show that the Wilson loop average in Yang-Mills theory is rewritten 
into a bosonic string theory. We see that 

J d 4 x(e, upa d»B^) 2 = J d A x2{-B, v d 2 B pv - B^d^B^ + B up d,d x B Xu ) (2.51) 

using the integration by parts. There are two ways to perform the integration over 
B^v field to obtain the Wilson loop average. One way is to identify the theory Sd[B] 
written in terms of B with a gauge-fixed version of a gauge-invariant master theory 
with an action Sm[B,A], see Appendix B. This way is applicable to the case up to 
£ = 0. In this case, the gauge fixing condition for B can be chosen such that d p B pa = 
which greatly simplifies the calculation, see Appendix [B] (or Appendix D of |61|| ). 
Another way is to deal directly with the action Sa[B] without imposing any condition 
on B. This way is applicable to any i. In this case, the Wilson loop average has extra 
contributions corresponding to the boundary terms in the coordinate representation. 
However, it turns out that the boundary term is not responsible for the area decay, 
see Appendix 0. Therefore, we can neglect the boundary term to obtain the area 
decay of the Wilson loop average. 

Thus, the relevant part of the action responsible for the area decay of the Wilson 
loop average reads 



d 4 x - B I 

ARM 2 v vpa 



M 



,2\2 



M 2V 



B pa + 



(2.52) 



where 



M 



4/to//t 



M 



2k- 



K (0) 



a vc + 



KM 



«:= (1/k q + 1/kI)-\ (2.53) 



15 



In the Landau gauge, W^^^ M „ is proportional to the transverse projection 



d^d 2 <9„, see @H. Therefore, we can regard that the tensor structure of KjfJ is simply propor- 
tional to <5„„, since antisymmetry of h^ v implies dpd a h^ a = 0. 
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The general form of the full gluon propagator in momentum space is given by 
&£?(p) =^ B [(W 2 - P,PM(P) + >T l Pi jp v B{j>) + 5, v C(p)}- 1 (2.54) 



--5 



AB 



- P T + - P L 

p 2 A{p) + C{p) p 2 B{p) + XC{p) ^ 



[2.55) 



where P// u := 5^ v — p^Pu/p 2 and P^ u := p^p u /p 2 ■ The asymptotic form of the full 
propagator in the large p region reads in the Landau gauge 



V(p) = * 

while in the small p region 

^ B (p) = * 



1 Cjp) 
p 2 A(p) p 4 A 2 (p) 

1 p 2 A(p) 



0(1// 



P, 



0(p 4 



P T . 



(2.56) 



(2.57) 



The OPE result is inconsistent with the absence of both current and dynamical 
gluon mass, i.e., C(p) = 0. In fact, the OPE calculation with the renormalization 
group improvement (|2.38| ) in the Landau gauge A = yields 



3?(p) =K b (p) + w [ ^*{p)(\k ■ K (x)) 



(2.56 



1 / lnp/Ap 
\\n/i/A 



13 Nc 
' 8 Po 



+ 



N c g 2 (l^ 2 ) /lnp/Ac 



3 N c 

4 00 



2(N 2 ~1) p± 



ln/x/A 0/ 



+ 0(1/ p 6 

(2.59) 



P T . 



From the consistency of the OPE result with the general asymptotic behavior (|2.56j ), 
we are lead to the existence of non-vanishing gluon mass function C with the loga- 
rithmic behavior: 



13 N c 
6 Po 



C(p) 



1) K 2 ^ln/i/Ao, 



61 N c 
12 (3 



2(N 2 



(2.60) 



On the other hand, we find 

d A p 



~-a Q 



(2tt 
d A p 



4 o AB y-P 



)W^ A J(p)( 
G AB (-P)[^± B (P)-D 



AB 



(2.61) 



By substituting the asymptotic form fl2.57|) into ( 2.61|) , we obtain in the Landau gauge 
A = 



a,. 



~-a Q 



d A p 
N c g 2 



G 



AB, 



-p)8 



AB 



2{Nl 



1V2 



-i 



1 p 2 A(p) 
d A p 



+ 0(p 4 



(2tt)< 



G AB (-p)5 



AB 



'lnp/Ap' 
ln/i/Ao 



# 



0{p 2 



(2.62) 
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where G AB {p) is the Fourier transform of G AB (x — y) defined by ( p.32j ) and # is 
a numerical number to be calculated. The integration over four momenta in (|2.62|) 



gives a dimensionless number. Therefore, K^a vc is proportional to the inverse of the 
vacuum condensate^ 



K (o) 



Or. 



N c g 2 



2{Nl 



-Alt 
IV 2 



-\ -1 



2 ) 



(2.63) 



In the absence of gluon pair condensation (&Z 2 ) = 0, M defined by (|2.53|) vanishes, 
M — 0. Roughly speaking, M 2 is proportional to the gluon condensation of mass 
dimension 2, i.e., apart from a numerical coefficient 



M 



2 

g a vc . 



^2.64) 



In this paper we have taken into account only the first two terms in Q2.41 ) corre- 
sponding to £ = and £ = 1. An advantage of including the £ = 1 term in addition 
to £ = one is that the ultraviolet (UV) cutoff A can be removed in the final stage 
of calculating the string tension to obtain a finite value for the string tension (See 
Appendix 0). 

The integration over the field is achieved by Gaussian integration to yield 



W{C) ^exp 



J2 

^ i=l,2 



i) m (e,/( 



-5 2 + M^r 1 e) + ^(^) 2 (e,/e) 



where we have defined Q 



Ml 



L ~ / / M 2 



/ M 2 

X :=M\MH{Ml - M\) = M 2 /Jl - 4^. 



(2.65) 

(2.68) 
(2.69) 



We mention two special cases. We consider the first limit (apparent London limit): 
M < M -> oo, namely, M 2 = M -> oo and M 1 = M < oo. In this case, the mode 



16 Of course, it is possible to perform the loop calculation to obtain the coefficient. The RG 
improved OPE ( 2.60} ) leads to # = as shown above. However, in order to obtain reliable 



results over the whole momentum region, it is desirable to solve the flow equation of the non- 
perturbative renormalization group. The relevant results from this viewpoint will be reported in 
forthcoming papers |101|. 

17 We have defined the L 2 inner product: 



(Q,IKQ) := J d±x J d i y e fw (x)I fwpr7 K(x 1 y)e pa (y)- 



(2.66) 



In particular, 



(6,76):= / dSe^aOWQpcrO")- 



(2.67) 
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with a huge mass M 2 decouples from the theory and we have only to include the mode 
with a mass Mi alone corresponding to i = 0. On the other hand, in the second limit 
(apparent Bogomol'nyi limit) defined by M x = M 2 , namely, M x = M 2 = M/V2 ^ 
V2M, we must include the second term too, i.e., £ = 0, 1. 



c 



A X " 



T 



R 



Figure 2: A planar Wilson loop with side lengths R and T. 



Now we proceed to calculate the Wilson loop average. We choose a rectangular 
loop C with side lengths R and T (Fig. 0) for simplifying the calculation. It is shown 
by explicitly calculating the integral in the argument of the exponential ( |2.65| ) (see 
Appendix |B| or Appendix F of |JT|) that the Wilson loop average exhibits the area 
law decay : 

W(C) S exp[-^Area(C)], a st = J 2 g 2 k£- In ^, (2.70) 

where Area(C) = RT is the minimal area of the surface spanned by the rectangular 
loop C and a st is the string tension. The string tension in the first limit reads 
a st = J 2 g 2 k^- In ||, while in the second limit a st = J 2 g 2 M 2 . What the area law 
holds irrespective of the shape of the loop C can be understood through the string 
representation given below. 



2.8 Step 8: String representation 

Now we obtain the string representation of the Wilson loop average in Yang-Mills 
theory. Let o = (a 1 , a 2 ) be a two-dimensional coordinate on the world sheet, while 
the target space coordinates is expressed by = x M (a). Then the infinitesimal 
surface element, 

dS, lu {x{a)) = ^ g(a)tf, u (a)d 2 a, (2.71) 
is expressed by the so-called extrinsic curvature tensor of the surface, 

= 1 d a x^(a)d b x„(a), (2.72) 
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c 




dS(x) 
x= x (a) 




dS (x) 
x = x(a) 



Figure 3: Interacting two surface elements dS pu (x) and dS pa (x') in the surface E 
spanned by the Wilson loop C. 



where the determinant g(o~) = det ||5W>(o")|| is calculated from the induced met- 
ric tensor of the surface defined by g a b(cr) = d a x^{a)di,x^{a) with the derivative, 
d a = gf^ (a = 1,2). We use the conformal gauge for the induced metric gab(o~) = 

We can observe in the expression (|2.65|) that there are massive propagating modes 



mediating the interaction between two surface elements dS pu (x) and dS pa {x') in the 
surface E spanned by the Wilson loop C, see Fig. [5| It is shown that the derivative 
expansion in powers of Q := (a' — o~) a Mi leads to 

(6,/(-<9 2 + Mfy l Q) 

d 2 S^(x(a)) I d 2 S pr7 (x'(a'))I^ pa (-di 2 + Mf)-\x{a),x'{a')) 



Jd 2 a^g 



am; 



AO) 



±M-^} 2) g ab (d a t^)(d b t 



hi/ j 



+ 0(C 4 



(2.73) 



where is the n-th moment defined by 



(n) 



<? 1//4 C a , from the Green function G 



'-d 2 + M 2 )-h 



Jd 2 z(z 2 ) n G Mi (z 2 ) for zt :-- 



47T 2 \x— X 



™K x {M\x- 



x'\), with Ki being the modified Bessel function, see Appendix ||. 

Thus we obtain the confining string theory for the Wilson loop average: 



W(C)^exp(-S cs [x}), 

S cs =0St / d 2 o-y/g + ccq 1 / d 2 a y /gg ab d a t^ u d b t flu + K t d 2 o^fg~R + ■ 



(2.74) 



where[^] the first term is the Nambu-Goto action, the second term is the rigidity term, 
and the third term is the intrinsic curvature term and the parameters are determined 

18 Here we have fixed a surface spanned by the Wilson loop. Taking into account the summation 
over the possible surfaces, we obtain W[C) = J T>x l _ L (a)J[x] exp(— S cs [x}) with J[x] being the Jaco- 
bian associated with the change of variables from the vorticity tensor to the target space coordinates 
(see e.g., tiftM). 
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as 



47T 



<T 



-1 




(2.75a) 
(2.75b) 
(2.75c) 



" < 0, 

2 1 
Kt =-J 2 g 2 k — > 0. 

3 Air 

The details of the calculation is given in Appendix [E] (see also Appendix G of [|6l|j ) 
which could be compared with a review article [^, [F5[] where a similar calculations 
were performed. The string theory just derived is nothing but the rigid string with a 
negative rigidity term proposed by Polyakov |50| . (See and p9|- p2| , |4^ , ^9| , pl| ,[52 



for the detailed studies of the rigid string.) In the absence of the vacuum condensate 
of mass dimension 2, i.e., a vc — > 0, the string tension vanishes a s t —>■ and the area 
law decay is lost. In other words, non-zero vacuum condensation a vc ^ is suggested 
from the non-zero string tension a st ^ 0. 



2.9 Step 9: Dual Ginzburg-Landau theory 

We proceed to show that the gluon condensate of mass dimension 2 is equivalent to 
the monopole condensation which is believed to be a mechanism for the dual super- 
conductivity of gluodynamics.0 The dual transformation (or Fourier transformation) 
can be used to transform the massive antisymmetric tensor field theory into another 
equivalent theory. It is well known |76|,[T^1 that in D-dimensional space-time, a mass- 



less antisymmetric tensor field of rank p (p-form) is dual to the (D — p — 2)-form, 
while a massive antisymmetric field of rank p is dual to the (D — p — l)-form. In 
fact, it is shown in Appendix [FJ that the effective theory ( |2.48| ) retained up to the 
quadratic terms in the massive Kalb-Ramond [^] field B^ v is dual to the massive 
vector field theory described by the dual Ginzburg-Landau (DGL) theory in the 
London limit: 



Sdgl = I d x 



(2.78) 



where U^ u := d^U u — d u ll^, [/£, is the Dirac string tensor defined by U^ u (x) := 
JgQ^x) and the radial mode of the scalar field is frozen and the phase is absorbed 
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In the previous paper JlCj ] , we stressed the importance of the term: 

J d 4 xh, u f^ = J (PxB^fV = J d 4 xd,b„*f^ = J d 4 xb,k^, (2.76) 

where we have defined the monopole current: k^ — d v * fY v after the Hodge decomposition of = 
d l Ja v — d„b^ + ■ ■ ■ . The VEV of the quadratic term yields the mass term of the dual (Abelian) gauge 
field bf,, if {{k^{x)k v {y))v)YM = g 2 5^5 4 {x - y)K H as 

(((7 d^xh^jAj ) y ) YM = J d 4 x J d 4 yb tl {x)b v {y){{k^)ku{v))v)YM - J d^mgt/Jfc). 

(2.77) 

where 77i| is proportional to K. This fact leads to an interpretation that the dual gauge field becomes 
massive due to the dual Meissner effect caused by monopole condensation {{k^(x)k tl (x))y)YM 7^ 0. 
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into the vector field Up to make the vector field massive with the mass 



mu = k 1/2 M ^ \jg 2 ^ 2 ). (2.79) 

Here the kinetic term of corresponds to the mass term of and vice versa. 

In the limit M — > 0, the dual gauge mass vanishes mu — > 0. In this limit, if 
the field is rescaled as — > U^/mu, then the configuration U^ u = becomes 
dominant, so that the vector theory written in terms of = <9 M </> reduces to a scalar 
field theory written in terms of <fi. Therefore, for the dual Ginzburg-Landau theory 
to describe the massive dual vector as a consequence of the dual Meissner effect, the 
dynamical generation of the kinetic term of B^ v is indispensable through non-zero 
value of a vc . In other words, the mass gap caused by the gluon (and ghost) pair 
condensation in the original Yang-Mills theory is equivalent to the dual Meissner 
effect caused by monopole condensation in the dual theory. 

As has been shown, the London limit (an extreme limit of the type II dual su- 
perconductor) of the DGL theory is obtained by including the quadratic term in the 
effective theory. In order to reproduce the DGL theory beyond the London limit, it 
is necessary to include the quartic term in h, see Appendix |]. 

The final result still depends on k. It should be remarked that the original theory is 
independent of k. Therefore, if we can take into account all the terms of the cumulant 
expansion, the result should be ^-independent. However, the dual theory obtained 
from the truncated effective action may depend on k. We hope that the dependence 
will be substantially small. If so, the truncation of the infinite series turns out to be 
a good approximation. The precise estimation of this dependence will be given in a 



forthcoming paper based on the exact renormalization group method ||101| 



3 Remarks 

We can consider an alternative step for Step 5. This enables us to discuss the relation- 
ship between the gauge potential correlator srf — srf and the field strength correlator 
T — T briefly. The effective theory Sd[h] is also obtained in connection with T — T 
correlation function. In fact, jY. v has also the manifestly gauge- and BRST-invariant 
form: 

fU x ) ■= V ( x ) • ^M - 9~ l V{x) ■ (DpVix) x D v V(x)), (3.1) 

where := — igl^, ■)■ It is easy to see that the first term V(x) ■ T^ix) and 
the second term f2 lu/ (x) := V(x) • (D fM V(x) x D v V(x)) = —Q v ^{x) of RHS of this 
equation is separately gauge invariant. Here note that (J~^ 1/ (x)J-' B (y))YM alone is not 
gauge invariant. Therefore, using this expression in ( |2.25|) , we obtain an alternative 
BRST-invariant action: 

s h =J d 4 *i- KM ~^J dix J d 'yK^)hAy)(UU x )fUv))T)TM + ■■■ 
= J d A x^-h^(x)h^(x) - ^g- 2 J d 4 x J dSKMKM(^,M^ P M)T) c YM 

- ^ J d 4 xj d'y{V\x)V B {y))T(^)TfM)YM KM KM + ■ ■ ■ ■ 

(3.2) 
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By expanding h(x) and h(y) in the formal power series in r := x — y around the 
coordinate X := (x + y)/2, the above action is cast into 



S h 



d 4 X 



#(i) 



+ 



(3.3) 



where we have assumed the translational invariance for the correlation function as 
functions ofr 2 = (x — y) 2 and we have defined 



■= / d\x-y) {{x-yfY(UUx)f V P M)T) 



v con 
'YM- 



(3.4) 



This effective action ( |3.2|) just derived should be compared with the previous one 
( 2.30|) . An advantage of this approach is that the equation preserves the manifest 
gauge and BRST invariance. The disadvantages are as follows. The evaluation of the 
field-strength correlator (e.g., (^t(^)^ r ^-(l/))yM) m the effective action ( |3.2|) is much 
harder than the gluon field correlator (^(i)^ B (|/))y^ in ( |2.30p . Furthermore, it 
is difficult to relate this version of the effective theory with the gluon condensation 
of dimension 4, i.e., (Fm{ x )J~£,(x))ym, since the OPE of the composite operator 



Fuv( x )3~m{y) does n °t include the gluon condensation of mass dimension 4 and hence 
cannot give the manifestly gauge and BRST invariant result, see [71 1. However, the 
following identification of the gauge-invariant correlation function might be related 
to the approach of the stochastic vacuum model (SVM) .PI 



((C(x)fV(y))v 



com con 
/YM 



rDi{{x - yf) + [d^x - y) p 5 ua - (x - y) 
+ d u ((x - y) a 5 w - (x - y)p6f^)]D 2 ((x - y) 2 



cr w upj 



(3.5) 



Once the two functions D\ and D 2 are determined, it is possible to calculate Kffi 
of ( |3.4p to fix the action of the effective theory (p.3[). 

In the Fock-Schwinger (FS) gauge x^srf^x) = (or coordinate gauge) [82] - 



it 



is possible to compare both approaches, since the gauge field ^ is written in terms 
of the field strength T^ v in the FS gauge: 



x 





oo 



da ax v T^ax) 



-E 



n=0 \ n + 2 ) n 

1 



x v x^x^ ■ ■■x"«[D lJl (0), [AJO), [■ ■ • [Ajo),:?>] •■■]]] 

1 



- o^(0)^ - -iVPa(0),^(0)] - -x u x x x T [D x (0), [D T (0),F^(0)] + ■ • • 
16 o 

(3.6) 



20 In the non-perturbative study of QCD, the cumulant expansion is extensively utilized by the 
SVM [^H where the different version of the non-Abelian Stokes theorem is adopted. In the SVM, 
the approximation of neglecting higher order cumulants is called the bilocal approximation. The 
Gaussian correlator has been calculated on a lattice [f79[-|3lf and the validity of bilocal approximation 
in SVM was confirmed by Monte Carlo simulation on a lattice The relationship between two 
gluon correlators and the dual field correlator in the supposed DAH model has been investigated 
in |80| . The authors would like to thank Dmitri Antonov and Nora Brambilla for this information. 
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where -D M (0) := — zg[^(0), •]. This implies the gluon propagator 



1 

4 X y (AT2 _ \)D(D - 1) 1 Ml/ pCT + 



(3.7) 



Hence, two types of gluon condensations are related as 

K^x) • ^(x)) = -^(FapiO) ■ ^(0))x 2 + • • • . (3.8) 

It should be remarked that the FS gauge does not retain the translational invariance. 
In fact, we do not assume translational invariance from the beginning in this paper 
The translational invariance implies that the VEV of O is equal to the VEV of the 
integrand Q(x): 

(O) = Q- 1 f d 4 x(Q(x)) = Q- 1 f c/ 4 x(Q(0)> = (Q(0)). (3.9) 



In the translationally non- invariant gauge, the situation is quite different. In the 
FS gauge, the vacuum condensate of mass dimension 2 is related to that of mass 
dimension 4 as 



(O) = ir 1 J d A x (^^{x) ■ s/^x)} = g^(^(o) ■ ^(o))^- 1 J d 



4 x x 2 + ■ 



(3.10) 



In order to make this quantity well-defined, it is necessary to introduce the long- 
distance or IR cutoff .0 

In the Lorentz gauge, defining 



0(q) = Q- 1 J d A x e iqx 



-^(x)^(x) + \iV(x)V{x) 



(3.11) 



we find that the vacuum condensate of mass dimension 2 proposed in JL2}] is equivalent 
to the statement for condensation of the zero-momentum mode (O) = (0(q = 0)) ^ 0. 
Therefore, this can be interpreted as the simultaneous Bose-Einstein condensation 
of gluon pair and ghost-antighost pair. The phenomenological description of gluon 
pair condensation has been discussed in |87|-|89|. This issue is to be investigated 
subsequently in detail. 



21 The authors in |8(| have concluded that the condensate (g/^ ) is zero, provided that the trans- 
lational invariance holds and the approximation of the vacuum dominance in intermediate state is 
accepted. Therefore, our results suggest that this approximation is not so good. Contrary to the 
claim in |^6| , the vanishing of the condensate is not essential to retain the gauge invariance of the 
vacuum energy. 

22 If the theory is defined on a finite space-time volume, then (O) is automatically finite. However, 
it does not guarantee the existence of the infinite volume limit. 
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4 Conclusion and discussion 



In this paper we have proposed a strategy to obtain the effective field theory or the 
effective string theory of the Yang-Mills theory with an insertion of the Wilson loop 
operator. In fact, we have obtained a number of effective theories which are equivalent 
to the original Yang-Mills theory within the validity of the approximations adopted 
here (weak- field approximation and lower derivative expansion). 

The first effective theory is written in terms of an antisymmetric tensor field h 
which is coupled with the surface (vorticity tensor) spanned by the Wilson loop C. 
At first, the antisymmetric tensor h was introduced as an auxiliary field without its 
kinetic term. After integrating out all the other fundamental fields which are expected 
to become massive, h has acquired the kinetic term due to the gluon pair condensation 
in the Landau gauge (a special case of the simultaneous Bose-Einstein condensation 
of gluon pairs and ghost- ant ighost pairs in the general Lorentz gauge). This effect has 
been estimated as a non-vanishing vacuum condensate of a composite operator with 
mass dimension 2 which is a color-singlet quantity and a BRST invariant combination 
of the gluon field and the ghost (antighost) field [|13j. 

The second effective theory is a confining string theory with a rigidity term, which 
was derived through the non-local action for the vorticity current tensor. The bosonic 
string theory known as a rigid string was conjectured to be an effective theory of 
gluodynamics a long time ago. At least, in the weak-field approximation and the 
derivative expansion, our result confirms this conjecture.^ 

The third effective theory is the dual Ginzburg-Landau theory obtained through 
the duality transformation in the sense of the electric-magnetic duality. The mass 
gap in the dual theory is understood as a manifestation of the dual Meissner effect 
caused by the monopole condensation. 

An advantage of our approach is that we must assume from the beginning neither 
the existence of mass gap nor a bare mass term in Yang-Mills theory. Such assump- 
tions were often adopted to study the non-perturbative feature of Yang-Mills theory. 
Just as the dynamical generation of quark mass associated with the spontaneous 
breaking of chiral symmetry, the gluon and ghost masses are of dynamical origin 
(namely, generated as a consequence of the dynamics of the relevant theory from con- 
sistency) associated with the breaking of the scale symmetry which is however broken 
anomalously. 

Another advantage is that the vacuum condensate a vc of mass dimension 2 in 
Yang-Mills theory is calculable analytically p^, TE, 72 and also numerically |T9|.|20 



The non- vanishing vacuum condensate stems from the gluon pairing |M],|91] at least 
in the Landau gauge. In the gauges other than the Landau, we need to take into 
account the ghost-antighost pairing [Tj] , |T5| ] to maintain the BRST invariance, as 



pointed out in |Tj| . It is also interesting to examine the confinement in quantum 



electrodynamics (QED) due to photon pairing which is expected to occur in 

the strong coupling phase 



Thus the dual Ginzburg-Landau theory is able to describe the dual superconduc- 



23 In the modified MA gauge, all order contributions from Vl^ v were incorporated by making use 
of the dimensional reduction to the two-dimensional nonlinear sigma model and instanton calculus 
in the dilute gas approximation within the framework of a deformation of a topological field theory 
as a reformulation of Yang-Mills theory |23] . 
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tivity of the Yang-Mills vacuum. The type of dual superconductor can be determined 
by specifying the action of h more accurately. The quadratic action in h corresponds 
to the London limit, an extreme case of the type II, where the radial degrees of free- 
dom of the scalar field is freezed and only the phase degrees of freedom is active. 
Within the quadratic approximation, we have obtained a finite value of the string 
tension after removing the ultraviolet cutoff. The removal of the cutoff was made 
possible first in this paper by taking into account the next-to-leading order of the 
expansion which was overlooked so far. 

To go beyond the London limit, it is necessary to specify the quartic term which 
is also calculable from the Yang-Mills theory as a higher-order cumulant according to 
the strategy indicated in this paper. In order to perform the calculation, however, we 
need to know more detailed properties of the correlation function of the gluon field 
(and ghost field). 

The integration of massive modes or high-energy modes can be interpreted as 
a step of Wilsonian RG. For a justification of this viewpoint, however, we need to 
perform more serious studies of the Wilsonian RG, as has been tried in [p5| , |99| , p!00 |. 
This issue will be investigated also in a subsequent paper [|101|| . 

Moreover, it is not yet clear whether the large N (color) limit simplifies this 
project and is able to give more clear connection between the Yang-Mills theory and 
the string theory. This issue is to be investigated in future. 
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A Cumulant expansions 

The cumulant expansion up to the order Kq 2 reads 



m ^]=i(/A(I/ v /X-i/ M ;)>7 

+ 7^2 ( / dA% J d4 y (2 Kv( x )flv{x) - ^fvv( X )f»v( X )) 

/I 1 \ \ con 
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The integration by parts allows us to rewrite the space-time integral of hn U fX u into 



^x\h^(x)f^{x) = I d i x[d u h^ u (x)}\V{x) ■ ^{x)\ 



2 



-~ I d*xg- 1 h li ,,l,r)V(.r)-(j) !l V(.r) xO„\'(.r)). (A.2) 



First, we consider the contribution to the double cumulant expansion from m = 1: 
g{F[h, £/})ym- In what follows, we omit all the terms which do not couple to 
field. Hence, the first term in ( |A.1| ) reads 



com con 



d 4 x[d u h^(x)](V(x)) c v ™ ■ {^{x))Tm 



1 ' d 4 xg-% v (x)(V{x) ■ (d^V(x) x d v V(x)))™ 



Kq 

= -^rj rf 4 ^- 1 V(x)(^ 1/ (x)) c v ? n , (A.3) 

where we have used (V(x))y n = (^(x))ym = and defined 

%u(x) := V(x) ■ (d^V(x) x d v Y{x)). (A.4) 
Similarly, the second term in ( A.l ) has the expectation value: 



d'x f d*yd„h, u (x)dph a p(y)(V A (x)V B (y))™(^(x)^ (y))™ 



+ 1T29~ 2 I d ' x I d 4 yh^x)h a p(y)(n, u (x)n a p(y)) c °^ (A.5) 

OKq J J 

where we have used (^WlyM = again. 

Next, we consider the contribution from m = 2: ({F[h, £^}) 2 )ym- By taking into 
account ( |A.3| ), the contribution comes from 



((/ rf4 ^V( x )/^( x )} ) 

d 4 x J d 4 yg'% u (x)h a p(y){n^{x))^(n a p(y))^ (A.6) 
We can put (^^(x))^? 11 = from Lorentz invariance. Thus we arrive at ( |2.30|) . 



B Manifestly covariant quantization of antisym- 
metric tensor field 

In this appendix, we adopt the Minkowski formulation. 
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B.l Massless case 



We discuss the gauge fixing of a second-rank antisymmetric tensor gauge field A^ v 
whose Lagrangian in Minkowski space-time is given by 



£ = -\{e^d v A"' 



■a\2 



This Lagrangian is invariant under the hypergauge transformation, 

SA^x) = d^ v {x) - d v ^(x). 
In order to fix the gauge, we adopt the gauge fixing condition for A^, 



d u A 



0. 



(b.i; 



(B.2) 



(B-3) 



Then the gauge fixing (GF) and Faddeev-Popov (FP) ghost term is obtained based 
on the prescription of Kugo and Uehara |)3| as 



-id* 



C^d^A^ + ^-B u ) 



(B.4) 



where we have introduced the vector ghost C, antighost C and the Nakanishi-Lautrup 
(NL) field B. Note that C and C are independent fields and that C* = 0,0^ = C . 
Then the nil potent BRST transformations are defined by 



d^C v (x) 



d B A^(x) = 

S B d(x) =0, 
6 B B^(x) =0. 



d u C^(x), 



Hence, the explicit form of C\ reads 



c x = B'd^v + icr^c - d v c»\ + y 



(B.5) 



(B.6) 



However, the Lagrangian C\ and hence Co + C\ is still invariant under the trans- 
formation of the vector ghosts C M and C^, i.e., 5C M (x) = id^9(x),5C^(x) = id^^x). 
Therefore we must fix the gauge for the vector ghost and antighost. Here we consider 



the gauge fixing conditions, <9 M C M 
GF+FP term, 



and d^C,, = 0. Thus we obtain an additional 



Co 



46 



13 



d{d"C^ + a 2 P) 



iS 



13 



N (>C M + a 3 B^ 



(B.7) 



where the nil potent BRST transformations of the additional fields are supplemented 

as 

S B N(x) =P(x), 
6 B P(x) =0, 
6 B d{x) =5 (1) (x), 
6 B B (1) (x) =0. (B.8) 
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field 


rank 


ghost number 


A 


2 





C 


1 


1 


d 





2 


C 


1 


-1 


B 


1 





N 








P 





1 


d 





-2 







-1 


A 


1 





a 





1 


a 





-1 


B' 









Table 1: The ghost number of the field with the indicated rank. 



The explicit form of C 2 reads 

C 2 = -iB^d^C^ - ia A B {l) P + dd^d^d - iPd^C^ + Nd^B^ (B.9) 

where we have defined := a 2 — 03. Note that P and B^ anti-commute. For the 
assignment of the ghost number of each field, see Table. 1. Two vector fields C M and 
are primary ghosts, and three scalar fields d, d and N are secondary ghosts. Three 
fields B^ P and are the Lagrange multiplier fields for the condition, d v A^ = 0, 
8^0^ = and d^C^ = 0, respectively. Thus we obtain the GF+FP term for the 
Lagrangian C , i.e., Cgf+fp = Ci + C 2 . Now all the gauge degrees of freedom are 
fixed. Thus the full Lagrangian density is given by 



£-tot — £0 + £-gf+fp, 

- iB^d^Cf, - ia 4 B w P + dd^d^d - iPd^C^ + Nd»B^. (B.10) 



The massless antisymmetric tensor field stands for the massless spin-0 field as a 
physical mode. It is possible to show that all the unphysical modes decouple leaving 
correctly one physical mode ||102|. PI 



24 The above result is the summary of the results obtained by Townsend [102], Kimura [102] and 
Hata, Kugo and Ohta [104]. The same result can be obtained within the framework of the extended 
theory for the constrained system based on the canonical Hamiltonian formalism on the extended 
phase space, the so-called the Batalin-Fradkin-Vilkovisky (BFV) formalism, see e.g. the original 
papers and a review [ 105 1 . 
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B.2 Massive case 

Next, we consider the theory of an antisymmetric tensor field with the mass term. 
The Lagrangian in Minkowski space-time is given by 

C%[A\ = -l^F*") 2 ~ \m 2 (A^. (B.ll) 

This Lagrangian with the mass term is no longer invariant under the hypergauge 
transformation of A^ v . However, the invariance is recovered by introducing an addi- 
tional vector field A M in such a way that 

C?[A, A] = -\{e^d v A^f - \{mA, v + d,A v - d y A,f. (B.12) 

Actually, this Lagrangian is invariant under the combined transformation, 

5A^ v (x) = d^ v (x) - d v ^(x), dA^x) = -m^(x). (B.13) 

Moreover, it has another invariance under the transformation, 

SA^(x) = 0, 5A^x) = d»u{x). (B.14) 

Therefore, we define the BRST transformation of A and A as 

SbA^x) =d li C v (x) - duC^x), (B.15) 

S B A^x) = - mC^x) + d^C'ix), (B.16) 

where we have introduced the ghost C and an extra ghost C The BRST transfor- 
mations of C and C are determined by the nil potency as 

SbC^x) =idpd{x), (B.17) 

d B C'(x) =imd(x), (B.18) 

S B d(x) =0. (B.19) 

The BRST transformation of antighosts C and C' is defined by 

SbC^x) =^(x), (B.20) 

SbB^x) =0, (B.21) 

5 B C'(x) =iB'{x), (B.22) 

S B B'(x) =0, (B.23) 

6 B d(x) =5 (1) (x), (B.24) 

6 B B {1 \x) =0. (B.25) 

Moreover, we introduce N and P in such a way that 

6 B N(x) =P(x), (B.26) 

d B P(x) =0. (B.27) 



28 



Thus the nil potent BRST transformation is determined for all the fields. In order to 
fix the gauge degrees of freedom, we must add the GF+PF term Cqf+fp- 
A good choice isQ 



C 



GF+FP 



l8j 



a i 



C v d'Ap, - d v N - aA v + -^B v + d (d"C„ + bC + a 2 P) 



a 



(B.28) 



fJbV 



d u N - aA v + -j-B v 



+ iC^id^Cu - d v C^) - d v P - a{d v C - mC v )} 
- iB^^C^ + bC + a 2 P) + dd^d^d + bmdd 



a: 



+ B\d»A, + -B<) + iC'd^d.C - mCp), 



(B.29) 



where ai,a 2 ,a' are gauge fixing parameters and a, b are parameters specified later. 
Roughly speaking, this corresponds to the gauge fixing condition, d^A^ = 0, d^C^ = 
= d^C^ and M A^ = 0. For the total Lagrangian C' tot := Cq' + Cqf+fp, the 
generating functional of the theory is given by 

Z := J VA^VA^VB^VC^VC^VdVdVNVPVB^VC'VC'VB'exp i J d 4 xC' tot . 

(B.30) 

After integrating over B and B', the sector containing A and A reads 



with 



Z := J VA^VA^VNexp i J d 4 x£' tot 



C'U =C^[A] - l -mA, v {d,A v - d v Afj) - \{d,A v - d v A»? 



(B.31) 



_ —id^AA 2 - — 
2a' K m; 2ax 



(B.32) 



since other fields decouple from the relevant sector. Then, performing the integration 



over N, we obtain 



Z=J VA^exp {i J d 4 x (C™[A] - J- (d»A^f 
x J PA^expjz J d 4 x£™[A,A] 



fB.33) 



where we have defined 

[A, A] := - KiA^K - d„A„) - ](d,A u - d u A 



4 



a , . . 2 a 



y (A,y + -(d»A, u )A u + — 
2a\ oti 2a>i 



(d^AJA-^d'Av). 



(B.34) 



25 The authors would like to thank Atsushi Nakamura [ 106 1 for helpful discussions on this Ap- 
pendix. 
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If we choose a = ma\, the cross term {d v A tlv )A tl cancels with mA llu (d f j,A u — d v A^). 
Hence, C™[A, A] becomes independent of A field: 



£™[A) = --m 2 (d lt A v -d u A ll ) 



aim 



(Am) 2 + 



aim 



(d»AJA- L (d»A u ). (B.35) 



Thus the A field decouples from the theory of A. Consequently, the theory with 
C' tot [A,A] := Cq' + £gf+fp is equivalent to the theory with C X [A\ := C$[A] - 
2^- l&'A/u,) 2 . In particular, the choice of the Landau gauge a\ — yields 



Z = J VA.J^A^exp {i J d A xC Q n [A] 



^B.36) 



On the other hand, if we choose a = 0, the theory with C' tot [A, A] := C^' + Cqf+fp 
is equivalent to the theory with C%\A\ := C™'\A, A] — (<9"yl M j,) 2 . In particular, the 
choice of the Landau gauge a x = yields 



Z = J DA p PA^(^ F )exp ji | ^HAA] 



^B.37) 



Note that the massive antisymmetric tensor gauge theory stands for the massive 
spin-1 theory [|76[ 107] . The antisymmetric tensor A^ v = —A Vil has six components 
in four-dimensional space-time. The gauge fixing condition d v A^ = yields three 
independent relations among the six components, since a trivial condition d tl d u A IJjV = 
is satisfied for the antisymmetric tensor A^ u . Therefore, A^ v has three independent 
components. This number agrees with the degrees of freedom for the massive vector 
field of spin one. 



Moreover, it is possible to consider a simpler gauge ||1 OS 

Cgf+fp = - iS B [C^A^ + d& 

=B fl A fl + iC^dyP' - mC^) - %B {l) C + mdd. 



(B.38) 
(B.39) 



This corresponds to the gauge fixing condition, A^ = and C = 0. The integration 
over N, P and B' is trivial, since CJ tot does not include them. The B integration leads 
to the constraint 5(A M ), 



Z = J VA^VA^A^) J VC^VC^VdVdVB^VC'VC' exp i J d A x£" c 
where 

£' = £™>[A, A] + i&iPyp - mC^) - iB {l) C + mdd. 



(B.40) 



(B.41) 



When we consider the sector of A and A, the sector described by other fields decouples 
and we obtain 



Z = J VA^VA^AJ expjz J A£HAA]} = / PA^expjz J d 4 xC^[A]j . 

(B.42) 

Therefore, we recover the original theory which is written by the A field only with 
the Lagrangian C™. 
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C Integration over antisymmetric tensor field and 
the boundary term 

By including the source term, the Lagrangian is given by 
£>d[B] =—B ilv B ilv + — (—B^d^B^ - B^ u d u d x B Xfl + B l/ll d ll dxBx v ) 



1 „ _ , 1 



— -^(—B^uB^ — B^du—dxBx^ + B Ufi d^—d\B\ u ) — -^-^B^I^^d B p(T 

— [ — B ^ u (d 2 ) 2 B ^ — B^ u d 2 d u dxBxfj, + B up d 2 d tJi d\Bx u ] + i-B^J^, 

(C.l) 

where J^ v := gJ *Q^ U . The saddle point equation of the Lagrangian in the momentum 
representation is given by 

-B^u{p) + (\- -o-o) \p 2 Bhv{p)+PvP\B\h{p) -PnP\B\v{p)\ 



p 2 _ 1 

+ —IuvpaBpa(p) + —p 2 \p 2 B^{p) + PuP\B Xfl (p) - PnP\B Xu (p)] + iJpAp) = 0. 



(C.2) 



Here we introduce the projection operators I and P defined by 



Pfiua/3 ■ ^{T^T^p T^pT ua ) , . 5fj, u 2 , 

z p 



where 

\ivp<jl paa/3 ^pvafii PpupaPpaafS Ppuafii v^"^) 

1 pvpcrPpcrctfl Ppupa^ puafi Ppvctfl- (^"^v 

Note that the projection operators have the properties: 

IpuafS Iup.a/3 lafSpu Ipvftai i.^^) 

Ppvafl Pup.a/3 PafSpu Ppvfia-, (^-'•8) 

and 

P Pl Inv P ivP P2I a(3p 2 a =^{PpPa9ufi - PpPfSQua - PuPa9p.fi + PuPfSQpa) 

=^(/-PW (C.9) 
For an arbitrary antisymmetric tensor A^ u , it follows that 

IpvapA^ = A^ u , A^ V I^vaf} = A a/3 , (C.10) 

p 2 Ppva(3A aP = p 2 A /IU - PvPpAtf + PpPliA v p. (C.ll) 
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Thus the saddle point equation is rewritten as 
1 p 2 



1 J ) I^vctP + 



1 . £. . t \ P 



(C.12) 



It is shown using the properties of the projection operators that the inverse of the 
operator [al + bP] is (a + 6) _1 [/ + (b/a)(I — P)\. Hence the saddle point value of B is 



kM\M 2 



{p 2 + Mf)(p 2 + Ml 



1 + M' 2 



(I-P) 



fiuaft 



-iJapip)), 

(C.13) 



where we have defined M 2 := rf / M 2 := 7 4 / ? ? 2 > M' 2 := c<j 2 /k, k := + 1/k 2 ) 1 , 
and (7 4 = ~kM 2 M 2 ) 



(M. 



2,1, 



2 y U 2 



1 ± 



\ 



1 - *7- 4 ( 4> + "4 



(C.14) 



By substituting this value back into the Lagrangian in the momentum representation, 
the B integration is performed to obtain 



ex PS - T 



d A p . 



ip(x-y) 



4 J (2tt) 4 



J^(-p)t a p(p) 



X 



X 



v 2 + Ml p 2 + Mf 



=1 + pL(-\ _i_ pf^ 

1 T" M ,2 



(I-P) 



(C.15) 



If we neglect the term proportional to (I — P), this agrees with the argument of the 
exponential of ( p. 65 ). It turns out that the same result is obtained by the action ( |2.52 ) 
which is obtained from ( U.lj ) by setting d\B\^ = 0. The London limit is obtained by 
taking the limit 7, u — > 00. 

It is possible to show that the term proportional to (/ — P) yields the boundary 
term. The identity ( |C9| ) implies 



exp< 



X 



X 



p 2 + Ml p 2 + Ml 



=1 + vL (1 + vL 

k ~ M 2 V M 2 
1 4- 



exp< 



(I — P)nua(3 

d 4 p 



X 



X 



(2tt) 4 \p 2 + Ml p 2 + Ml 



+ 2g v p&°d y a H(x,y) 



where 

H(x,y) :-- 



X 



X 



d 4 p I 

(2tt) 4 \p 2 + Ml ~ p 2 + Ml) f' 



-1 



^ M 2 



(l + ^ 



1 4- -2— 

1 ^ M' 2 



M 2 J c ip(x-y) 



(C.16) 



(C.17) 
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Hence the second term in the exponential is written using the partial integration into 
the boundary term: 

~ J d^xdlJ^x) J d A yd v j: au {y)H{x,y) = — ^g 2 J 2 dx u dy u H(x, y). (C.18) 

D Calculation of the Wilson loop average 

We choose a rectangular loop with side lengths R and T in the X\ — x 4 plane: 



e^(z) = 6 IIl 5^5(z 2 )5(z 3 )e(z l )e(R - Zl )e(z 4 )e(T - * 4 ). 

Then the Fourier transformation is obtained as 

e^(p) = J d 4 ze^(z) e - ip - z 

S^5 u4 [ R dz ie - ipiZl [ T dz 4 e- ip4Z4 

-e 2 sin : 



(D.l) 



Jo Jo 

, , 2 -viR . piR 2 -PiT _ p 4 T 
Onid^—e 2 sm 



Pi 2 p A 

In the momentum representation, we have 

d A p 



(D.2) 



(9, I(-d 2 + Mf)^Q) = J W(-^ 2 + Mf)-\p)Q aP {-p). (D.3) 

For large L = R and L = T, we can apply to (p.2| ) the formula: 

lim f— V = 7rL5(a), (D.4) 

Then ( p.3j ) reads 



(e,/(-a 2 + Mf)- 1 e)-/ ? ^ I (2vr) 2 ™(p 1 )^4)(-9 2 + M2)- 1 (p) 



(27T) 



(27T) 2 

d 2 p 



-# + M 2 )- 1 (0,p 2 ,p 3 ,0) 
1 



i = 1,2 . 



(D.5) 



(27r)> 2 +p| + M 2 

Each integral is logarithmically divergent. By taking into account two terms, the 
logarithmic divergence of the integral is removed after introducing an ultraviolet cutoff 
A as 

E (-l) m (©> H-d 2 + M 2 )- 1 ®) =TR I - ^ 

i=l,2 



(27T) 



1 



=Ti? lim 



pi+pi + Mf pi + pi + Mi 
d\ P f 



A^oo Jo 47T 



H 2 + M 2 \p\ 2 + M| 



!rfl ili m ln A2 + M ? A « 



4ttA^ - Mf A 2 + M| 

_ 1 , M 2 
=Ti?— In— -. 
2vr Mi 



(D.6) 
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Similarly, we can calculate the inner product: 

(0,10) := J d 4 xO^(x)I^ pa O pa (x) = 
d A p 



d 4 p 



(27T) 



-(27TyTR5( Pl )5(p 4 ) 



--TR 



d 2 P _ r „A 2 
(2tt) 2 4tt' 



(D.7) 



This part corresponds to the quadratic divergence (i.e., the volume divergence in 
two-dimensional space). This should be removed by the renormalization prescription. 

Therefore we obtain the A-independent and finite result. Thus the Wilson loop 
average exhibits the area law decay, 



T 2 2~_X M 2 



W(C) S exp(~a st RT), a st = J g R-^ In 



(D.8) 



E Derivative expansion of the non-local string ac- 
tion 



In this section, we begins with the expression 

1 



W(C) =exp 



—J 2 g 2 k [ dS^(x) [ dS^(y)G IMUpa (x,y) 



G 



fiupcr 



%y) :=l 



X 



X 



-A + Mf -A + Mi 



where x '■= MiMf/(Mf — M 2 ). We define the Euclidean propagator: 



G M (x] 



:-a e + m 



2\-l, 



x,0) 



d 4 k 



ik-x 



The propagator is obtained in the closed form: 



G M {x) 



1 M 



At: 2 



x\ 



MMlx]) 



(2vr) 4 k 2 + M 2 ' 
• 2 ^i(M/0 



4tt 2 



-M z 



(E.l) 
(E.2) 

(E.3) 
(E.4) 



where Ki(z) is the modified Bessel function and £ is the correlation length defined 
by £ = M _1 . In fact, ( E.4|) is obtained as follows. Substituting the identity, 

1 '^dse~ s ^ +M2 \ (E.5) 



k 2 + M 2 jo 

into ( E.3|) and performing the Gaussian integration over the four momenta k, we 
obtain 



Gm(x) = / dse 
Jo 



dse 



-sM 2 



-sM 2 



d A k 



-sk +ik-x 



exp 



x 

Is 



(27T) 



167T 2 Jo 



c?s-r exp 



-sM z 

4s 



(E.6) 
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The above result (|E.4) is immediately obtained by applying the integration formula 



I09fl : 



/ ds s u ~ l exp 
Jo 



-^-M 2 s 
s 



2((3/M 2 ) vI2 K u (2^3m' 2 ) (3fy3 > 0, KM 2 > 0)), 

(E.7) 



to the case v — — 1 and /3 = x 2 /4, since K_ v {z) = K u (z). 
In Euclidean space, 



G^upoix, X ) — I^upcrX [(*Mi (x — x) — Gm 2 (x — X )] . 



We define 



7,:=/ dS^(x(a)) I ^(xK))WG Ml ((x(a)-xK)) 2 ^ 



It is shown |26|,|T5f that the derivative expansion of Ji in powers of 
leads to 

1 



(E.8) 
(E.9) 

(E.10) 
(E.ll) 



Ji = J d 2 a^g [4C 2 ^ (0) - -^^ 2) 9 a \d a t, v ){d b t, v ) 
where no summation is understood over i and we have defined the moment, 

Jtt ] := I d 2 z(z 2 TG Mt (z 2 ), zf := gW$. (E.12) 



Here we have used the conformal gauge for the induced metric, g a b(&) = \j g(c)^ab 

which leads to ( a ( b g a b = g~ 1 ^ 2 g a bZ a z b = z a z b 5 a b '■= z 2 . Thus, the confining string 
theory derived in this paper is characterized by the parameters, 



ost =J 2 g 2 ~KX J d 2 z[M^ 2 G Ml (z 2 ) - M^ 2 G M2 (z% 

"o 1 = " I^V^X / d 2 zz 2 [M^G Ml {z 2 ) - M^G M2 {z 2 )], 



Kt =lf A J 9 *X / d 2 zz 2 [M^G Ml (z 2 ) - M^G M2 (z)\. 



(E.13a) 
(E.13b) 
(E.13c) 



By substituting (|E.4j) into ([E.13a ), we obtain 



=J 2 g 2 k^ 



2it\z\d\z 



2it\z\d\z 



Ki(\z\) 



;e.i4) 



where we have used K\{x) = —K' (x) and K (oo) = and introduced an ultraviolet 
cutoff A (since K Q (0) = oo). Note that the asymptotic behavior of the modified Bessel 
function Kq(z) for z <C 1 is given by 



z 2e~ lE 
K {z) = -{ lE + \n-)=\n , 



(E.15) 
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with 7s being Euler's constant 7^ = 0.5772 ■ • • . Thus, for sufficiently large A, we 
obtain the A- independent finite result. Incidentally, the asymptotics of K p (z) for 
z > 0, 



K p (z) 



;e.i6) 



means that K p (z) decreases exponentially for large z. After removing the cutoff A, 
the above expression reduces to a finite value, 



. (M 2 
2?r a VMi/ 



(E.17) 



This is an advantage of our approach over the previous one [p7| . The coefficient of 
the rigidity term is calculated as 



On 



\z 


2 m 


A) 


Ml 


z 





\z\d\z\ 



\z 


2 K, 


(1 


A) 




z 





2„2~ X 



- J 9 K 



647T 2 



An An 

mJ~mJ 



j2 2~ 

- J g k 



16n 



< 0, 



where we have used the integration formula (v = l,fJ>= 3), 



dxx^ l K v [ax) = 2^ 2 a~ / T 



\L — V 



r 



fi + V 



(9fy/ > »!/). 



(E.18) 



(E.19) 



Here note that the integral in ( |E.18| ) is finite and we don't have to introduce the 
cutoff. Similarly, the parameter n t is calculated as 



K t = -o a o = o J 9 > 0. 

3 3 107T 



(E.20) 



F Path-integral duality transformations 

F.l Dual Ginzburg-Landau theory in the London limit 

We introduce the field strength H of an antisymmetric tensor field B (the so-called 
Kalb-Ramond field) by H := dB, i.e., 



H, 



d x B^ + d u B uX + d u B 



(F.l) 



Note that the field strength H is invariant, H^x — > H^x, under the hypergauge 
transformation, B — > B + d(, i.e., 



5 



B^ ■= B 



(F.2) 



26 The string tension is a free energy per unit length of the string. It is well known that the free 
energy has a logarithmic dependence in the Ginzburg-Landau theory. The London limit corresponds 
to Mi — > 00. Hence, the string tension reduces to the expression, a s t = J 2 g 2 H^Ko ("x") , since 
K (M 2 /A) •(). 5 
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We require the invariance of the measure T>B^ V under the hypergauge transformation. 
However, the Lagrangian Cd[B] does not have the invariance under the hypergauge 
transformation due to the existence of the mass term (B^) 2 . Nevertheless, we can 
recover the hypergauge invariance of the antisymmetric tensor theory by introducing 
a new vector fieldQ A M which transforms as 

A, - A£ := A, - 0, (F.3) 

In fact, the combination B A := B + dA, i.e., 

B A V = B^ + - d v A^ (F.4) 

is invariant under the combined transformations, ( |F.2p and ( |F.3j ). Therefore, the 
Lagrangian, 

C m [B, A] := C m [B A ] = + M A„ - d v Arf + J-^„ A ) 2 + J-^tf^) 2 , 

V 7 (F.5) 

is also invariant, i.e., C m [B,A] = C m [B^ , A^}. 

As we have recovered the hypergauge invariance, we need to fix the hypergauge 
invariance in quantizing the dual magnetic theory C m [B,A]. For this purpose, we 
adopt the gauge-fixing condition, 

d v B^ = 0, (F.6) 

for the antisymmetric tensor field.0 Under this condition, the derivative terms of 
B^ u in C m [B A ] reproduce the corresponding terms of £d[B], 

{H^f = - W^8 2 B^ - Gd.B^B^ - -3B^ u d x d x B^, 
(d x H Xfll/ ) 2 ={d X d x B, u + d,d x B uX + d u d x B x J 2 - B, u (d x d x ) 2 B^. (F.7) 

Therefore, C m [B,A] under the condition (|F.6|) reduces to 



C m [B, A] = + 0„A„ - d u AJ 2 - -^B^d x d x B^ + -LB^(d x d x ) 2 B^ . 

V 7 (F.8) 

It is shown that the theory with C m [B, A] reduces to the original theory given by Cd[B] 
by integrating out A field after fixing the gauge freedom of A M , see Appendix [B]. 

Thus we obtain an alternative dual description of low-energy Gluodynamics in 
terms of B^ v and A M . Especially, for G = SU(2), we obtain 

W(C) = Zll J VB^B^) J VA„ exp {-S M [B A ; C}} , (F.9) 



27 This field plays the similar role to the Stiickelberg scalar field in the massive vector theory which 
recovers the gauge invariance of the vector field. 

28 In the manifestly covariant quantization of the gauge theory, we need to introduce the ghost as 
is well known. However, it is not enough for the antisymmetric tensor gauge theory, since we need to 
introduce the ghost for ghost in order to completely fix the gauge degrees of freedom. Such a theory 
is called a reducible theory. In this subsection we treat the theory in a naive manner. However, the 
result is unchanged if we take into account the reducibility of the theory. See Appendix |b|. 
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where 



S M [B A ;C] = J d 4 xC m [B A ] + l -Jg jf dS^B, u . 

Apart from the third term in £ m [B,A], the above action QF.5|) coincides 
action of confining string proposed by Polyakov [50| in the weak field limit. 
We define 



S M [B A ; 0] = J d 4 x {c m [B A } + l -JgB, u *Q 

Then the Wilson loop average is given by 

W(C) = Z M [B A ;Q}/Z M [B A ;0}, 



(F.10) 
with the 

(F.ll) 

(F.12) 



where 



Z M [B A ;Q] 

VB^iVB^) J V(,exp{-S M [B<;e]}. 



12rf 



= J VB^J^B^expl- J d A x 

x exp j— i J d^x-JgB^y *0 A 
x Jv(,e W {-Jd*xj-(Bp 
Note that the path integral transformation holds, 
|l?^exp{-|d 4 x^(^ 



(H Xllu ) 2 + —(d x H 



X/ii/ ) 



= J Vi^e^d^-B^exp^- J d 4 xj-(£, ul ( . 
since the constraint, 

e^ pa d p (e, u - BfJ) = 0, 

is solved by 

BpjJ 3 U (^^y 2.6, 

Moreover, we introduce the auxiliary (Abelian) vector field by 



I = B^ 



where U^ v is the (dual) field strength defined by 



(F.13) 



;f.i4) 



(F.15) 



(F.16) 



5(e^d p (£^ - BfJ)) = J PC/ M exp l-i J d A x *t^(V " B ^)} > (F.17) 



(F.18) 
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By using the identity ( |F.17| ), the integration over in (|F.14|) can be performed as 



J VU^expli J d\ *U^B^ J m^expl-J 



d x 



_(7 ) 2 



VUf, exp 



d x 



k{U^) 2 - i *U IMU B i 



Hi/ 



(F.19) 



Hence the equality holds,^] 



PC, exp ^- / d 4 x-(BlY 



VUn exp 



d x 



(F.23) 



Thus, the theory is rewritten in terms of and B^ v as 

z M [u,B-e] 



x exp 



d x 



x exp j-i y d 4 x^JgB^ u *Q fll/ \ ■ 
By change of variable —* £7^ + ^JgQ^ u , we arrive at the expression, 

= | VU^exp i- J d A xK (u^ + l -JgQ, 
x JvB^S^B^expl- J d 4 x 



(F.24) 



1 



1 



12l , 2 -(^) 2 + ^p»*. 



(F.25) 



29 Another way of deriving the equality is as follows. The argument of the exponential in the LHS 



IS 



(B + d(,B + dQ = (B, B) + (B, d() + (dCi 5 ) + d C) ~ B) + {d(, dC), 



(F.20) 



under the condition 5i? = 0. The last term decouples after the Gaussian integration of (. On the 
other hand, the argument of the exponential in the RHS is cast into 



d*x [(CV) 2 - *U^{B^)} = (dU,dU) - (*dU,B) = (U,SdU) - (U,*dB). 



(F.21) 



Suppose the Lorentz type gauge condition SU = 0. We introduce the NL (zero-form) field <f>. Then 
the Gaussian integration over J7 M field yields 



(U, AU) - (U, *dB) - (SU 7 <f)) = (U, AU) - (U, *dB + #) 

Sd , 



1 Sd 
(*dB + dcf>, — * dB + #) = [B, —B) 



(B,B) 



(F.22) 



under the condition SB — 0. In this derivation, we must insert the constraint S(d^ i U tl ) in the measure 
VU^. The identity implies that there are many ways of extracting the transverse modes of B. 
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We change the variable B, v into the new variable Y~ M asQ 



| M^(d^)exp j-| d 4 x 
= J VB,J(d u B^) exp |- J d 4 x -ie^U,d u B pu + J-^(^) 2 + -l-^H^f 
= J VY^d^) exp |- J d'x -i2U,Y» + ^Y^ + -^(0„n - d v Y,f j , 

(F.26) 

since the constraint <9 M Y~ M = can be solved by an antisymmetric tensor field in the 
form, 



(F.27) 



The massive antisymmetric tensor field B^ denotes the massive spin-1 field Y~ M whose 
canonical mass dimension is three. This should be compared with the massless an- 
tisymmetric tensor field which stands for the massless spin-0 field, see [ 102 - 104 ]. In 
this step, the number of independent degrees of freedom is conserved, since Y~ M and 
B pv have three independent components. The path-integral duality transformation 
indicates the correspondence: 



Y f 



1 



c iiupa 



d v B 



1 



e^ pu E vpu = 3 V *B^ d v W v . 



2 o 

Furthermore, the integration over Y p is performed after introducing the new vari- 
able 9 to remove the delta function of the constraint d^Y^ - = 0, 

! "2 , 1 /o v a v ^2 



j V6 exp < 


-J d 4 x 


j V6 exp < 


- J d 4 x 



iBdJT* - 2iWU ll + — 2 YJ + —{d,Y v - d v Y, 

1 1 



T>8exp 



d 4 x 



4 



7 



r -~p ^ ^_ 5 2 +7 4/^2 



(F.29) 



Finally, after rescaling of — > -L^, we obtain the dual Abelian gauge theory, 

z M [u,e-,e] 



■■ J VU^l/p) I V6exJ - J d l 



X 



4 



(U^ + u 



S \2 
fiu) 



7 



+ 5 (^ + ^) _ y+7W i- 



/IV 



r 



(F.30) 



30 From SB^ = 0, there exists a three-form such that = 5Z& = <5 * WW = *eWW 
Then ff( 3 ) := dB^ = d*dW w = *5dW {1) = *rW, or = *ff( 3 ). Therefore, SY^ = 5*H& = 
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where is the Dirac string tensor, 

U%{x) :=J 9 e, u (x), (F.31) 
d v U%(x) =J* 4 := J g J^dr d -^5\x - x(t)). (F.32) 

Here we have inserted the delta function 5(8^11,) for fixing the gauge for U^. This 
model has the dual U(l) symmetry, say the magnetic U(l) m symmetry, 

U tt ^U tt + d^, 9^9-V. (F.33) 

For large 7, the second term has the expansion: 

2 \ _1 / ^2 \ 



\n 2 {u, + d,9) (1 - ^d 2 \ Lr - ^pp ) (u v + dj) 



\v 2 u, ( 1 + ^d 2 + ■ ■ ■ J cp + ^{d.ef + o(d') 



=\r?U»U» + \^U,d 2 U» + \rf{d,9f + 0(9 4 ), (F.34) 

where we have used d^U, = 0. Thus we arrive at an effective theory of the Yang-Mills 
theory, 

Zapegt[U,6;0) = J VU^dMVO expl - J d 4 x£ K [U,6) J, (F.35) 
£k[U, 9] :=- A (l + £\ (U, u ) 2 + \rfU^ + \rf{d,9f. (F.36) 



After putting k — 1, this model ( |F.36| ) just obtained should be compared with the 



London limit A — ► 00 of the dual Abelian Higgs model or the dual Ginzburg-Landau 
theory with the Lagrangian, 



CdqlP, 0] = \{U, V + Up 2 + - i2g m U,)<P\ 2 + A(|0| 2 - v 2 ) 2 (F.37) 



where g m is the magnetic charge subject to the Dirac quantization condition, 

9m9 = 4tt. (F.38) 
The London limit is equivalent to putting \4>(x)\ = v = const., i.e., <p(x) = v exp[i2g m 9(x)}: 

Cdgl[U] = \(U, U + U* u ) 2 + + d,9){U, + 8,9). (F.39) 

The effective theory ( |F.36| ) is of the same form as the London limit of the DGL 
theory, except for the renormalization of the kinetic term of the dual gauge field. The 
dual gauge field becomes massive, whereas the 9 field remains massless. This is rea- 
sonable, since the field 9 corresponds to the Nambu-Goldstone (NG) mode associated 
with the spontaneous breakdown of the magnetic U(l) symmetry. 
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The dual U(l) symmetry is broken in the London limit. This corresponds to the 
infinitesimally thin flux tube connecting the quark and anti-quark. In the London 
limit, the Higgs mass = 2y/\v diverges, i.e., = 0001 m^ 1 = 0. This is the 
extreme case of the type II superconductor where > mjj. It turns out that the 
mass rtiu of the dual gauge field is given by 77, 

ri = mu = 2g m v = —v. (F.40) 
9 

F.2 Beyond the London limit 

We consider a small variation from the London limit. First, we introduce the fluctu- 
ation field <p(x) around the absolute minimum |0(x)| = v: (f)(x) = v + T(p(x), where 
r := 1/A <C 1. Therefore, for (f)(x) = \4>(x)\e %e ^ x \ the covariant derivative reads 

:= - 2ig m U^ = [rd„<p + i(v + Ttp)^ - 2g m U ll )}e w , (F.41) 

which leads to the kinetic term: 

i|ZV>| 2 =~r 2 (<9^) 2 + i(v + r V ) 2 (d,9 - 2g m U,) 2 

= \r 2 {d^f + ^(v 2 + 2rvip + rV) W - 2g m U ll ) 2 . (F.42) 

Moreover, the potential reads 

A(|0| - v 2 ) 2 =\{\<f>\ - vfM + v) 2 = \TV(2r) + r V ) 2 

=4tv V + 4vr V + rV 4 - (F.43) 

Hence the DGL Lagrangian up to the order (p 2 reads 



C DGL [U, 0] =-{U, v + U* v ) 2 + -v 2 (d,9 - 2g m U,) 2 



1, , I 

4 

T 

~2 



.2 



+ ~(d^) 2 + rv(d,6 - 2g m U,) 2 p + 4rv 2 <p 2 + 0(<p 3 ), (F.44) 



Integrating out the fluctuation field <p, we obtain 



Cdgl[U, 6} =-{U, u + U* u ) 2 + -v 2 (d,6 - 2g m U, 



v 2 



-(d^O - 2g m U fl ) 2 (-d 2 + m 2 )-\d^e - 2g m U^) 2 , (F.45) 

where m 2 := 8v 2 /r = 8Xv 2 ^> 1 and r) := 2g m v. 

We can determine the antisymmetric tensor field theory whose dual version agrees 
with ( F.45|) . From the correspondence ([F.28 ) and a relation flF.27 ), 



Ua «-> — 
M 2i 



Y v , (FAQ) 
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the effective theory must have an additional term 

^ J d 4 x(H x , u ) 2 (-d 2 +m 2 )-\H 



72r) : 
m 

m 

72^f 



X/iu ) 



1 - 



d x(H Xflu ] 
J d 4 x(H x ^) 2 (l + 



d 2 



d 2 d* 
— + — + ■ 



d A x(H Xflu ) 2 

\h x * 



(F.47) 



This term is quartic in h with the fourth derivative and higher derivatives. Therefore, 
the next order term in the cumulant expansion 



4! 



f j4 }_i fV _ f j4 }_fV fV~' 
J a X 2* l ti.vJ fui J a x ^IfMuIfiv 



con\ con 



V I YM 



(F.48) 



is expected to generate the term (|F.47|) . In the tree approximation for the correlation 
of V, it is easy to see that the next order term involves the piece: 



C 2 
24«g 



d*x I dM<(x)<(y))YM(< A (x)K B (y))YM 



x (d ■ h),(x)(d ■ h)„(x)(d ■ h) p (y)(d ■ h) a {y) 

C ~ 1 f d'x f A((<(^)<(y))yAf) 2 (9-/ i ) p (x)(9-/ i ) p (y)(9-/i) (T (x)(9-/ i ) (7 (y) 



(F.49) 



24/416 

-^M' { °'l d ' x{a - hnx)(d -" nx) 

+ J_ s (i) f d 4 X(d ■ h) 2 (X)d 2 (d ■ h) 2 (X) + 
16 J 

where we have defined (d ■ h)^ := d u h pu and the coefficient is defined by 

8® := / d 4 r((^ A (x)^ B (y)) YM ) 2 (r 2 Y (£ = 0, 1, • ■ ■ ). (F.50) 

Indeed, this expansion reproduces the desired terms (|F.47|) . 

The detailed analysis of higher order terms will be given in a forthcoming paper. 
Moreover, an argument supporting the above result will be given there from the 



viewpoint of the exact Wilsonian RG [55 



43 



References 



[1] Y. Nambu, Strings, monopoles, and gauge fields, Phys. Rev. D 10, 4262-4268 
(1974). 

G. 't Hooft, in: High Energy Physics, edited by A. Zichichi (Editorice Compos- 
itori, Bologna, 1975). 

S. Mandelstam, Vortices and quark confinement in non-abelian gauge theories, 
Phys. Report 23, 245-249 (1976). 

A.M. Polyakov, Compact gauge fields and the infrared catastrophe, Phys. Lett. 
B 59, 82-84 (1975). 

A.M. Polyakov, Quark confinement and topology of gauge theories, Nucl. Phys. 
B 120, 429-458 (1977). 

[2] G. 't Hooft, Topology of the gauge condition and new confinement phases in 
non-Abelian gauge theories, Nucl.Phys. B 190 [FS3], 455-478 (1981). 

[3] A. Kronfeld, M. Laursen, G. Schierholz and U.-J. Wiese, Monopole condensa- 
tion and color confinement, Phys. Lett. B 198, 516-520 (1987). 

[4] T. Suzuki and I. Yotsuyanagi, Possible evidence of abelian dominance in quark 
confinement, Phys. Rev. D 42, 4257-4260 (1990). 

[5] S. Hioki, S. Kitahara, S. Kiura, Y. Matsubara, O. Miyamura, S. Ohno, T. 
Suzuki, Abelian dominance in SU(2) color confinement, Phys. Lett. B 272, 326- 
332 (1991), Erratum-ibid. B281, 416 (1992). 

J.D. Stack, S.D. Nieman and R. Wensley, String tension from monopoles in 
SU(2) lattice gauge theory, Eep^at7MDgOTj 
(1994). 



Phys. Rev. D 50, 3399-3405 



H. Shiba and T. Suzuki, Monopoles and string tension in SU(2) QCD, [hep- 
lat/940401], Phys. Lett. B 333, 461-466 (1994). 

[6] O. Miyamura, Chiral symmetry breaking in gauge fields dominated by 
monopoles on SU(2) lattices, Phys. Lett. B 353, 91-95 (1995). 
S. Sasaki and O. Miyamura, Topolog ical aspect of Abelian projected SU(2) lat- 
tice gauge theory, ||hep-lat/9811029[] , Phys. Rev. D 59, 094507 (1999). 
S. Sasaki an d O. Miyamura, L attice study of UA(1) anomaly: the role of QCD- 
monopoles, ||hep-lat/ 98 10039] , Phys. Lett. B, 443, 331-337 (1998). 

[7] A. DiGiacomo, Monopole condensation and color confinement, |[hep-| 



lat/ 98020081 ; ||hep-th/9603029H . 

M.l. Polikarpov, Recent results on the abelian projection of lattice Gluodynam- 



ics, |hep-lat/9609020||. M.N. Chernodub and M.I. Polikarpov, Abelian projec- 



tions and monopoles, [ |hep-th/9710205 [ 
G.S. Bali, The mechanism of quark confinement, ||hep-ph/9809351 



[8] Z.F. Ezawa and A. Iwazaki, Abelian dominance and quark confinement in Yang- 
Mills theories, Phys. Rev. D 25, 2681-2689 (1982). 

[9] T. Suzuki, A Ginzburg-Landau type theory of quark confinement, Prog. Theor. 
Phys. 80, 929-934 (1988). 

S. Maedan and T. Suzuki, An infrared effective theory of quark confinement 
based on monopole condensation, Prog. Theor. Phys. 81, 229-240 (1989). 
T. Suzuki, Abelian confinement mechanism in QCD, Prog. Theor. Phys. 81, 
752-757 (1989). 

S. Maedan, Y. Matsubara and T. Suzuki, Abelian confinement mechanism and 
the QCD vacuum, Prog. Theor. Phys. 84, 130-141 (1990). 



44 



[10] K.-I. Kondo, Abelian-projected e ffective gauge th eory of QCD with asymptotic 
freedom and quark confinement, | |hep-th / 9 709 1 09j ] , Phys. Rev. D 57, 7467-7487 
(1998). 

K.-I. Kondo, | frep-th/ 9803053] , Prog. Theor. Phys. Supplement, No. 131, 243- 
255. 

[11] H. Ichie and H. Suganuma, Mono poles and gluon fields in QCD in the maximally 
abelian gauge, ||hep-lat/980805j , Nucl. Phys. B 574, 70-106 (2000). 



[12] C. Becchi, A. Rouet and R. Stora, Renormalization of the Abelian Higgs model, 
Commun. Math. Phys. 42, 127-162 (1975); Renormalization of gauge theories, 
Ann. Phys. 98, 287-321 (1976). 

I.V. Tyutin, Lebedev preprint, FIAN No. 39 (in Russian) (1975). 

[13] K.-I. Kondo, Vacuum co ndensate of mass dimension 2 as the origin of mass gap 
and quark confinement, ||hep-th/ 01 052991 , Phys. Lett. B 514, 335-345 (2001). 
K.-I. Kondo, Vacuum condensate of mass dimension 2 as the origin of mass gap 
and quark confinement (A review), ||hep-ph /0110013jj . 

[14] M. Schaden, Ma ss generation in continuum SU(2) gauge theory in covariant 
Abelian gauges, [|hep-th/9909011|] , 3rd revised version. 



[15] K.-I. Kondo and T. Shinohara, Ab elian dominance i n low-energy Gluodynamics 
due to dynamical mass generation, |[hep-th/ 00041 58|| , Phys. Lett. B 491, 263-274 
(2000). 

[16] Ph. Boucaud, A. Le Yaouanc, J. P. Leroy, J. Micheli, O. Pene and J. Rodriguez- 
Quinter o, Consistent OP E description of gluon two- and three-point Green func- 
tions?, [|hep-ph/0008043|| , Phys. Lett. B 493, 315-324 (2000). 



[17] F.V. Gubarev, L. S todolsky and V.I . Zakharov, On the significance of the vector 
potential squared, ||hep-th/001005l1 , Phys. Rev. Lett. 86, 2220-2222 (2001). 



INp- 



[18] F.V. Gubare v and V.I. Zakharov, Emerging phenomenology of (v4^m) 
ph/OOlTOfl , Phys. Lett. B 501, 28-36 (2001). 

R. Akhoury and V. I. Zakharov, On nonperturbative corrections to the potential 
for heavy quarks, | |hep-ph/971(M87| , Phys. Lett. B 438, 165-172 (1998). 

M.N. Chernodub, F.V. Gubarev, M.V. Polikarpov and V.I. Zakharov, Confine- 
ment and short distance physics, [|hep-ph/0003006|] , Phys. Lett. B 475, 303-310 
(2000). 



[19] 



Ph. Boucaud, A. Le Yaouanc, J. P. Leroy, J. Micheli, O. Pene and J. Rodriguez- 
Quintero, Testing the La ndau gauge opera tor product expansion on the lattice 
with a (A 2 ) condensate, [|hep-ph/0101302| , Phys. Rev. D 63, 114003 (2001). 
F. De Soto and J. Rodriguez-Quintero, Notes on the det ermination of the Lan- 
dau gauge OPE for the asymmetric three gluon vertex, ||hep-ph/0105063| . 



[20] Ph. Boucaud, J. P. Leroy, A. Le Yaouanc, J. Micheli, O. Pene, F. De Soto, A. 
Donini, H. Mouta rde, J. Rodriguez-Quintero, Instantons and (A 2 ) condensate, 
||hep-ph/0203TT9H . 

[21] G. Curci and R. Ferrari, On a class of lagrangian models for massive and mass- 
less Yang-Mills fields, Nuovo Cimento 32A, 151-168 (1976). 

L. Baulieu and J. Thierry-Mieg, The principle of BRS symmetry: An alternative 
approach to Yang- Mills theories, Nucl. Phys. B 197, 477-508 (1982). 

L. Baulieu, Perturbative gauge theories, Phys. Reports 129, 1-74 (1985). 



45 



R. Delbourgo and P.D. Jarvis, Extended BRS invariance and OSP(4/2) super- 
symmetry, J. Phys. A: Math. Gen. 15, 611-625 (1982). 

[22] K.-I. Kondo, Yang-Mills theory as a deformation of topological field theory, 
dimensional reduction and quark confinement, | [hep-th/ 980 1024| | , Phys. Rev. D 
58, 105019 (1998). 

[23] M.J. Lavelle and M. Schaden, Propagators and condensation in QCD, Phys. 
Lett. B 208, 297-302 (1988). 

[24] K.-I. Kondo, T. Murakami, T. Shinohara and T. Imai, Renormalizing a Becchi- 
Rouet-Stora- Tyutin-invariant composite operator of mass dimension 2 in Yang- 
Mills theory, ||hep-th/01 11251 , Phys. Rev. D 65, 085034 (2002). 

[25] M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, QCD and resonance physics, 
theoretical foundations, Nucl. Phys. B 147, 385-447 (1979). 
M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, QCD and resonance physics, 
applications, Nucl. Phys. B 147, 448-518 (1979). 

[26] D. Antonov, Str ing nature of confinement in (non-) Abelian gauge theories, 
||hep-th/ 990920911 , Surveys High Energ. Phys. 14, 265-355 (2000). 

[27] K.-I. Kondo, A confining string theory derived from QCD, CHIBA-EP-131, 
||hep-th/01 100041] . 



[28] J. Polchinski, Strings and QCD, [|hep-th/92l0045 



[29] G. 't Hooft, A planar diagram theory for strong interactions, Nucl. Phys. B 72, 
461-473 (1974). 

[30] M. Kalb and P. Ramond, Classical direct interstring action, Phys. Rev. D 9, 
2273-2284 (1974). 

[31] D. Forster, Dynamics of relativistic vortex lines and their relation to dual the- 
ory, Nucl. Phys. B 81, 84-92 (1974). 

J.L. Gervais and B. Sakita, Quantized relativistic string as a strong coupling 

limit of the Higgs model, Nucl. Phys. B 91, 301-316 (1975). 

A.M. Polyakov, Gauge fields as rings of glue, Nucl. Phys. B 164, 171-188 (1979). 

I.Ya. Aref'eva, Quantum contour field equtions, Phys. Lett. B93, 347-353 

(1980). 

[32] M. Liischer, K. Symanzik and P. Weisz, Anomalies of the free loop wave equa- 
tion in the WKB approximation, Nucl. Phys. B 173, 365-396 (1980). 

[33] O. Alvarez, Static potential in string models, Phys. Rev. D 24, 440-449 (1981). 

[34] J.F. Arvis, The exact qq potential in Nambu string theory, Phys. Lett. B 127, 
106-108 (1983). 

[35] P. Olesen, Strings and QCD, Phys. Lett. B 160, 144-148 (1985). 

[36] A.M. Polyakov, Fine structure of strings, Nucl. Phys. B 268, 406-412 (1986). 

[37] A. Polyakov, Gauge Fields and Strings (Harwood Academic Publishers, Chur, 
1987). 



46 



[38] H. Kleinert, The membrane properties of condensing strings, Phys. Lett. B 174, 
335-338 (1986). 

H. Kleinert, Spontaneous generation of string tension and quark potential, Phys. 
Rev. Lett. 58, 1915-1918 (1987). 

H. Kleinert, Thermal deconfmement transition for spontaneous strings, Phys. 
Lett. 189, 187-190 (1987). 

H. Kleinert, Dynamical generation of string tension and stiffness in strings and 
membranes, Phys. Lett. B211, 151-155 (1988). 

H. Kleinert, Exact temperature behavior of strings with extrinsic curvature 
stiffness in the limit of infinite dimension. Thermal deconfmement transition, 
Phys. Rev. D 40, 473-490 (1989). 

E. Braaten, R.D. Pisarski and S.-M. Tse, Static potential for smooth strings, 
Phys. Rev. Lett. 58, 93-96 (1987). 

P. Olesen and S.-K. Yang, Static potential in a string model with extrinsic 
curvatures, Nucl. Phys. B 283, 73-92 (1987). 

E. Braaten and C.K. Zachos, Instability of the static solution to the closed 
string with rigidity, Phys. Rev. D 35, 1512-1514 (1987). 

F. David and E. Guitter, Rigid random surfaces at large d, Nucl. Phys. B 295 
[FS21], 332-362 (1988). 

F. David, Random surfaces and the statistics of membranes, Phys. Reports 184, 
221-227 (1989). 

J. Polchinski and Z. Yang, High-temperature partition function of the rigid 
string, Phys. Rev. D 46, 3667-3669 (1992). 

M.I. Polikarpov, U.-J. Wiese and M.A. Zubkov, String representati on of 



the Abelian Higgs theory and Aharonov-Bohm effect on the lattice, [|hep- 
lat/9303007|l , Phys. Lett. B 309, 133-138 (1993). 



K. Lee, Dual formulation of cosmic strings and vortices, Phys. Rev. D 48, 2493- 
2498 (1993). 



P. Qrland, Extrinsic curvature dependence of Nielsen-Olesen strings, [[hep- 
th/940414q| , Nucl. Phys. B 428, 221-232 (1994). 



M. Sato and S. Y ahikozawa, Topological formulation of effective vortex strings, 
||hep-th/ 94062081 Nucl. Phys. B 436, 100-128 (1995). 

H. Kleinert and A.M. Chervyakov, Evidence for negative stiffness of QCD flux 
tubes in the large-N limit of SU(N), Phys. Lett. B381, 286-290 (1996). 
H. Kleinert, A.M. Ch ervyakov and V.V . Nesterenko, Evidence for negative stiff- 
ness of QCD strings, ||hep-th/9601030[ . 

A.M. Polyakov, Confining strings, ||hep-th/ 9607049] , Nucl. Phys. B 486, 23-33 
(1997). 

M.C. Diamantini, F. Quevedo and C.A. Trugenberger, Confining strings with 
topological term, ||hep-th/96l2T03H , Phys. Lett. B 396, 115-121 (1997). 

M.C. Diamantini and C. A. Trugenberger, Surfaces with long range correlations 
from noncritical strings, f|hep-tn/9712008|1 , Phys. Lett. B 421, 196-202 (1998). 
M.C. D iamantini and C.A. Trugenberger, Geometric aspects of confining 
strings, ||hep-th/9803046;1 , Nucl. Phys. B531, 151-167 (1998). 



47 



M.C. Diamantini, H. Kleinert an d C.A. Trugenbe rger, Strings with negative 
stiffness and hyperfine structure, fhep-th/9810f 71| , Phys. Rev. Lett. 82, 267- 
270 (1999). 

M.C. Diamantini, H. Kleinert and C.A. Trugenberger, Universality class of con- 
fining strings, ||hep-tli/99032US| , Phys. Lett. B 457, 87-93 (1999). 
M.C. Diamantini and C.A. Trugenberger, QCD-like behavior of high- 
temperature confining strings, | |hep-th / 0202 1 78| ] . 

M.C. Diamantini and C.A. Trugenberger, Confining strings at high tempera- 
ture, ||hep-th/ 0203051 . 

[53] A.M. Polyakov, String theory and quark confinement, | hep-th /971 1 002 1 , Nucl. 
Phys. (Proc. Suppl.) 68, 1-8 (1998). 

[54] E. Alvarez, C. G omez and T. Ortin, String representation of Wilson loops, 
||hep-th/ 980607511 , Nucl. Phys. B545, 217-232 (1999). 

[55] U. Ellwange r, Confinement, monopoles and Wilsonian effective action, |[hep- 
ph/9710326l Nucl. Phys. B 531, 593-612 (1998). 



U . Ellwanger, Field strength correlator and an infrared fixed point of the Wilso- 
nian exact renormalization group equations, Eur. Phys. J. C7, 673-683 (1999). 
U. Ellwanger, Monopole condensation and antisymmetric te nsor fields: Com - 
pact QED and the Wilsonian RG flow in Yang-Mills theories, [ |hep-th / 9906061|| , 
Nucl. Phys. B 560, 587-600 (1999). 

U. Ellwanger, M. Hirsch and A. Weber, Flow equations for the relevant part of 
the Yang-Mills action, Z. Phys. C 69, 687-697 (1996). 

[56] D. Antonov and D. E bert, String repre sentation of field correlators in the dual 
Abelian Higgs model, ||hep-th/ 9806151 , Eur. Phys. J. C 8, 343-351 (1999). 

[57] D. Antonov and D. E bert, Dual formu lation and confining properties of the 
SU(2)-Gluodynamics, [|iep-th/9902177|| . 

[58] R. Parthas arathy, Infrared region of QCD and confining strings, |[hcp- 
th/990325l . 



R. Parthasarathy and O. Suzuki, Abelianization of Wilson loop of non- Abelian 
gauge theory on complex manifold and Polyakov's conjecture, Preprint, June 
2001. 

V.I. Zakharov, Anatomy of a confining string, Physics Reports 320, 59-78 
(1999). 

[60] M. Baker and R. Steinke, A n effective string theory of Abrikosov-Nielsen-Olesen 
vortices, ||hep-ph/9905375| , Phys. Lett. B 474, 67-72 (2000). 
M. Ba ker and R. Steink e, Effective string theory of vortices and Regge trajec- 
tories, ||hep-ph/ 0006069| , Phys. Rev. D63, 094013 (2001). 

[61] K.-I. Kondo, Dual superconductivity, m onopole condens ation and confining 
string in low-energy Yang-Mills theory, |[hep-th/0009152|| (revised version in 
preparation). 

[62] K.-I. Kondo and Y. Taira, Non-Abelian Stokes Theorem and Quark Confine- 
ment in SU(3) Yang-Mills gauge theory, [ |hep-th/9906T29] , Mod. Phys. Lett. A 
15, 367-377 (2000); 

K.-I. Kondo and Y. Taira, Non-Abelian Stokes Theorem and Quark Confine- 
ment in SU(N) Yang-Mills gauge theory, |hep-th/ 991 1242|| , Prog. Theor. Phys. 
104, 1189-1265 (2000). 



48 



[63] T. Kugo and S. Uehara, General procedure of gauge fixing based on BRS in- 
variance principle, Nucl. Phys. B 197, 378-384 (1982). 

[64] D.I. Diakonov and V.Yu. Petrov, A formula for the Wilson loop, Phys. Lett. B 
224, 131-135 (1989). 

D. Diakono v and V. Petrov, Non-Abelian Stokes theorem and quark- monopole 
interaction, |hep-th/9606104j . 



[66 
[67 



[68 



[69 



[70 



M.B. Halpern, Field-strength and dual variable formulations of gauge theory, 
Phys. Rev. D 19, 517-530 (1979). 

N.E. Bralic, Exact computation of loop averages in two-dimensional Yang-Mills 
theory, Phys. Rev. D 22, 3090-3103 (1980). 

I.Ya. Aref'eva, Non-Abelian Stokes formula, Theor. Math. Phys. 43, 353-356 
(1980) (Teor.Mat.Fiz.43:lll-116 (1980)). 

Yu.A. Simonov, Cluster expansion, nonAbelian Stokes theorem and magnetic 
monopoles, Sov. J. Nucl. Phys. 50,134 (1989) (Yad. Fiz., 50, 213 (1989)). 
B. Broda, A four- dimensional Yang-Mills approach to polynomial invariants of 
links, Phys. Lett. B 262, 288-293 (1991). 

B. Broda, A topological-field-theory approach to the non-Abelian Stokes theo- 
rem, J. Math. Phys. 33, 1511 (1992). 

F.A. Lunev, Pure bosonic worldline path integral representation for fermionic 
deter minants, non-Abel ian Stokes theorem, and quasiclassical approximation in 
QCD, ||hep-th/9609166l Nucl. Phys. B494, 433-470 (1997). 



M. Hirayama and S. Matsubara, St okes theorem for loop variables of non- 
[|hep-th/9712120|| , Prog. Theor. Phys. 99, 691-706 (1998). 
Kanno, M. Ueno and H. Yamakoshi, Non-Abelian Stokes 



Abelian gauge field, 
M. Hirayama, M. 



frep-th/ 9806098| . 



theorem for loop variables associated with non-trivial loops, 
Prog. Theor. Phys. 100, 817-830 (1998). 

R.L. Karp, F. Mansouri, J.S. Rno, Product integral repres entations of Wil - 
son lines and Wilson loops, and non-Abelian Stokes theorem, [ |hep-th / 9903221|| ; 
R.L. Karp, F. M ansouri, J.S. Rno , Product integral formalism and non-Abelian 
Stokes theorem 



hep-th/99l0i75| , J.Math.Phys. 40, 6033-6043 (1999). 



K.-I. Kondo, Ab elian magnetic monopole dominance in quark confinement, 
||hep-th/980515l1 , Phys. Rev. D 58, 105016 (1998). 



M. Hirayama and M. Ueno, Non-Ab elian Stokes theo rem for Wilson loops as- 
sociated with general gauge groups, ||hep-th/9907063|| , Prog. Theor. Phys. 103, 
151-159 (2000). 

K.-I. Kondo, Quark confinement and deconfine ment in QCD fro m the viewpoint 
of Abelian-projected effective gauge theory, ||hep-th/9810167|| , Phys. Lett. B 
455, 251-258 (1999). 

K.-I. Kondo, A formulation of the Yang-Mills theory as deformation of a topo- 
logical fiel d theory based on the background field method and quark confinement 
problem, ||hep-th/9904045|1 , Intern. J. Mod. Phys. A 16, 1303-1346 (2001). 



D. Diakonov and V. 
lat/0008004 |. 



Petrov, On the non-Abelian Stokes theorem, Rhcp- 



[71] M. Lavelle and M. Oleszczuk, The operator product expansion of the QCD 
propagators, Mod. Phys. Lett. A 7, 3617-3630 (1992). 

[72] H. Verschelde, K. Knecht, K. Van Acoleyen and M. Vanderkelen, The non- 
perturbative groundstate of QCD and the local composite operator A^, [[hep- 
th/01050T8fl , Phys. Lett. B 516, 307-313 (2001). 



49 



[73] E.T. Akhmedov, M.N. Chernodub, M.I. Polikar pov and M.A. Zu bkov, Quantum 
theory of strings in an Abelian Higgs model, ||hep-t h / 9 50 50 70| ] , Phys. Rev. D 
53, 2087-2095 (1996). 

E. T. Akhmedov, Fermionic string from Abelian Higgs model with monopoles 
and 0-term, ||hep-th/96052T4]] . 

M.N. Chernodub, M.l. Polikarpov, A.I. Veselov and M.A. Zubkov, Strings and 
Aharonov-Bohm effect in Abelian Higgs model, ||hep-lat / 9804002)] , Phys. Lett. 
B 432, 182-186 (1998). 

[74] J. Polchinski and A. Strominger, Effective string theory, Phys. Rev. Lett. 67, 
1681-1684 (1991). 

[75] D.V. Antonov, D. Ebert and Yu.A. Simonov, Curvature expansion for the back- 
ground-induced gluodynamics string, | |hep-th / 9605086| j , Mod. Phys. Lett. A 11, 
1905-1918 (1996). 

D.V. Antonov and D. Ebert, Curvature expansi on for the gluody namics string 
including perturbative gluonic contributions, [[hep-th/9707097|1 , Mod. Phys. 
Lett. A 12, 2047-2056 (1997). 

[76] H. Kawai, A dual transformation of the Nielsen-Olesen model, Prog. Theor. 
Phys. 65, 351-364 (1981). 

F. Q uevedo and C.A. Trugenberger, Phases of antisymmetric tensor field theo- 
|hep-th/9604196H , Nucl. Phys. B 501, 143-172 (1997). 



nes, 



[77] P. Orland, Instantons and disorder in antisymmetric tensor gauge fields, Nucl. 
Phys. B 205 [FS5], 107-118 (1982). 

[78] H.G. Dosch, Gluon condensate and effective linear potential, Phys. Lett. B 190, 
177-181 (1987). 

H.G. Dosch and Yu.A. Simonov, The area law of the Wilson loop and vacuum 
field correlators, Phys. Lett. B 205, 339-344 (1988). 

[79] A. Di Giacomo, H. Panagopoulos, Field strength correlations in the QCD vac- 
uum, Phys. Lett. B285, 133-136 (1992). 

M. D'Elia, A. Di Giacomo and E. Meggiolaro, Field strength correlators in full 
QCD, [hep-lat/970503], Phys. Lett. B 408, 315-319 (1997). 
A. Di Giacomo, M. D'Elia, H. Panago poulos and E. Me ggiolaro, Gauge invari- 
ant field strength correlators in QCD, [|hep-lat /9808056H . 

[80] M. Baker, N. Brambilla, H.G. Dosc h and A. Vairo, F ield strength correlators 
and dual effective dynamics in QCD, ||hep-ph/ 9802271 , Phys. Rev. D 58, 034010 
(1998). 

[81] G.S. Bali, N. Bramb illa and A. Vairo , A lattice determination of QCD field 
strength correlators, [|hep-lat/9709079|1 , Phys. Lett. B 421, 265-272 (1998). 

[82] P. Pascual and R. Tarrach, QCD: Renormalization for the Practitioner, Lecture 
Notes in Physics Vol. 194 (Springer, Heidelberg, 1984). 

[83] M.A. Shifman, Wilson loop in vacuum fields, Nucl. Phys. B 173, 13-31 (1980). 

[84] V.A. Novikov, M.A. Shifman, A.I. Vainshtein, V.I. Zakharov, Calculations in 
external fields in quantum chromodynamics: technical review (abstract operator 
method, Foch-Schwinger gauge), Fortsch. Phys. 32 585- (1985) 

[85] V. Elias, T.G. Steele and M.D. Scadron, (qq) and higher- dimensional- 
condensate contributions to the nonperturbative quark mass, Phys. Rev. D 
38, 1584-1605 (1988). 



50 



I. Schmidt and J. -J. Yang, QCD condensate contributions to the effective quark 
potential in a covariant gauge, Phys. Rev. D 60, 114005 (1999). 

L.S. Celenza and CM. Shakin, Description of the gluon condensate, Phys. Rev. 
D 34, 1591-1600 (1986). 

L.S. Celenza and CM. Shakin, Phenomenological potentials for QCD parame- 
ters, Phys. Rev. D 35, 2843-2848 (1987). 

L.S. Celenza, C.-R. Ji and CM. Shakin, Coherent-state representation for the 
QCD ground state, Phys. Rev. D 36, 895-901 (1987). 

J.M. Cornwall and A. Soni, Couplings of low-lying glueballs to light quarks, 
gluons and hadrons, Phys. Rev. D 29, 1424-1434 (1984). 

H. Gomm, P. Jain, R. Johnson and J. Schechter, Scale anomaly and the scalars, 
Phys. Rev. D 33, 801-812 (1986). 

R. Fukuda, Tachyon bound state in Yang-Mills theory and instability of the 
vacuum, Phys. Lett. B 73, 33-38 (1978). 

R. Fukuda and T. Kugo, Dynamical theory of the Yang-Mills fields. I - 
Instability of the vacuum and singlet pair condensation-, Prog. Theor. Phys. 
60, 565-583 (1978). 

V.P. Gusynin and V.A. Miransky, On the vacuum rearrangement in massless 
chromodynamics, Phys. Lett. B 76, 585-588 (1978). 



K.-I. Kondo , Existence of confinement phase in quantum electrodynamics, |[hep- 
th/9803133|| , Phys. Rev. D 58, 085013 (1998). 

R. Fukuda, Photon pairing and the strong-coupling phase of massive quantum 
electrodynamics, Phys. Rev. Lett. 63, 482-485 (1989). 

T. Inagaki, M. Komachiya and R. Fukuda, Photon pairing and the strong cou- 
pling phase of massive quantum electrodynamics Cooper equation, Mod. Phys. 
Lett. A 5, 381-390 (1990). 

T. Inagaki and R. Fukuda, Photon pairing in massless quantum electrodynam- 
ics, Phys. Lett. B 273, 519-527 (1991). 

V.A. Miransky, Dynamics of spontaneous chiral symmetry breaking and the 
continuum limit in quantum electrodynamics, Nuovo Cimento 90A, 149-170 
(1985). 

W.A. Bardeen, C.N. Leung, S.T. Love, The dilaton and chiral symmetry break- 
ing, Phys. Rev. Lett. 56, 1230-1233 (1986). 

C.N. Leung, S.T. Love and W.A. Bardeen, Spontaneous symmetry breaking in 
scale invariant quantum electrodynamics, Nucl. Phys. B 273, 649-662 (1986). 

K.-I. Kondo, T. Mino and K. Yamawaki, Critical line and dilaton in scale- 
invariant QED, Phys. Rev. D 39, 2430-2433 (1989). 

T. Appelquist, M. Soldate, T. Takeuchi and L.C.R. Wijewardhana, in Pro- 
ceedings of the 12th Johns Hopkins Workshop on Current Problems in Particle 
Theory, edited by G. Domokos and S. Kovesi-Domokos (World Scientific, Sin- 
gapore, 1988). 

F. Frei re, Abelian proje cted SU(2) Yang-Mills action for renormalization group 



flows, ||hep-th/0110241|l , Phys. Lett. B 526, 405-412 (2002). 



51 



[100] U. Ellwanger an d N. Wschebor, Solvable confining gauge theories at large N, 
[[hep-tri/ 010709jJ|| , Phys. Lett. B517, 462-470 (2001). 

U. Ellwanger arid N. Wschebor, Confinement with Kalb-Ramond fields, | [hep- 
| th/OlOTT^ , JHEP 0110, 023 (2001). 

[101] T. Imai, H. Kato, K.-I. Kondo, T. Murakami and T. Shinohara, in preparation. 

[102] T. Kimura, Antisymmetric tensor gauge field in general covariant gauges, Prog. 
Theor. Phys. 64, 357-360 (1980). 

T. Kimura, Quantum theory of antisymmetric higher rank tensor gauge field in 
higher dimensional space-time, Prog. Theor. Phys. 65, 338-350 (1981). 

[103] P.K. Townsend, Covariant quantization of antisymmetric tensor gauge fields, 
Phys. Lett. B 88, 97-101 (1979). 

[104] H. Hata, T. Kugo and N. Ohta, Skew-symmetric tensor gauge field theory 
dynamically realized in the QCD U(l) model, Nucl. Phys. B 178, 527-544 (1981). 

[105] E.S. Fradkin and G.A. Vilkovisky, Quantization of relativistic systems with 
constraints, Phys. Lett. B 55, 224 (1975). 

LA. Batalin and G.A. Vilkovisky, Relativistic S matrix of dynamical systems 
with boson and fermion constraints, Phys. Lett. B 69, 309-312 (1977). 
Y. Igarashi, BRST quantization of the constrained system by Batalin-Fradkin- 
Vilkovisky, Soryushiron-Kenkyu (Kyoto) (in Japanese) 82, 414-436 (1991). 

[106] A. Nakamura, private communications. 

[107] A. Smailagic a nd E. Spallucci, The dual phase of massless/massive Kalb- 
Ramond fields, ||hep-th/0l06T73l , J.Phys.A: Math.Gen. 34, L435-L440 (2001). 

[108] C. Bizdadea and S.O. Saliu, The BRST quantization of massive abelian two- 
form gauge fields, Phys. Lett. B 368, 202-208 (1996). 

[109] I.S. Gradstheyn and I.M. Ryzhik, Table of Integrals, Series and Products (Aca- 
demic Press, Bosoton, 1980). 



52 



